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We  describe  a simple  expansion  of  the  Krylov-type  for  the  estimation  of  the 
electronic  correlation  energy  of  the  ground  state  of  closed-shell  systems  and  show 
its  relationship  with  Brillouin-Wigner  Perturbation  Theory.  We  demonstrate  that 
Brillouin-Wigncr  Perturbation  Theory  can  be  quickly  solved  in  a non-iterative 
fashion. 

Examples  are  presented  using  Intermediate  Neglect  of  the  Differential  Overlap 
Model  Hamiltonian  (INDO/1 ) which  indicates  that  this  scheme  captures  well  over 
95%  of  the  correlation  energy  captured  by  the  Configuration  Interaction  over 
double  excitations  method  (CID)  in  a fraction  of  the  computational  effort. 

We  also  described  a generalization  of  the  ground  state  method  to  obtain 
electronic  spectra  as  a single  excitation  on  top  of  the  ground  stale  Krylov-type 


ethod  (CIS)  expanded  wii 


average  over  doubles  and  triple  exciialions. 

Examples  are  presented  using  Intermediate  Neglect  of  the  Differential  Overlap 
Model  Hamiltonian  (INDO/1)  which  indicates  that  this  scheme  is  competitive  for 
transition  energies  and  moments  when  compared  with  the  much  more  demanding 
Configuration  Interaction  over  all  excitations  method  (CI-FuII). 

We  also  investigated  time  independent  response  of  a wavcfunction  constructed 
as  a Configuration  Interaction  over  singles  (CI-S)  using  the  first  order  Rayleigh- 

momcnt  and  polarizabilities  are  very  encouraging. 


CHAPTER  1 
INTRODUCTION 

Motivation 

Much  effort  has  been  and  is  being  expended  to  obtain  the  Potential  Energy 
Surfaces  (PES)  of  molecular  systems  as  a function  of  the  nuclear  coordinates. 
The  interest  in  the  calculation  of  Potential  Energy  Surfaces  is  rooted  in: 

i.  The  solution  of  the  Electronic  Hamiltonian  for  certain  nuclear  configurations, 
which  means  the  calculation  of  the  potential  generated  by  the  electrons  for 
the  nuclear  motion,  is  a well-defined  problem. 

ii.  It  is  a starting  point  for  the  understanding  of  a large  variety  of  problems,  c.g. 
structure,  dynamics  and  spectroscopy,  and  thus  of  all  chemistry. 

Here  we  will  concentrate  on  applications  to  electronic  spectroscopy,  photo- 
chemistry. and  time-independent  ground  state  properties.  So  we  are  interested  not 
only  in  good  transition  energies,  but  also  in  reliable  wave  functions  and  the  PES 
for  the  slates  of  interest. 

Section  3 of  this  chapter  presents  the  Electronic  Hamiltonian  we  are 


: attempting 


In  chapter  2 we  review  some  of  the  approximations  to  solve  the  Electronic 
Hamiltonian.  We  point  out  that  within  the  Spin-Symmetry  Adapted  Coupled 
Cluster  we  differ  from  Nakatsuji’s  work  1 1—1 8J  by  keeping  the  energy  evaluation 
consistent  through  second  order.  For  closed-shell  systems,  we  develop  diagram- 
matic rules  to  obtain  all  the  expressions  over  molecular  orbitals  instead  of  the 
usual  ones  over  spin  orbitals. 

In  chapter  3 we  analyze  a Krylov-like  expansion  for  a closed-shell  ground 
state.  We  show  that  this  expansion  is  closely  related  to  a perturbation-variational 
treatment  with  respect  to  Brillouin-Wigner  Perturbation  Theory  (BWPT).  This 
development  differs  from  other  perturbation-variational  methods  (e.g.  the  one 
by  Bartlett  and  Brhndas  [19-23])  as  it  does  not  use  an  expansion  over  positive 
powers  of  the  Hamiltonian.  We  develop  a scheme  that  uses  a power  series  in  the 
denominator  which  forces  a fast  convergence  of  the  perturbation  series.  Rather 
we  show  the  relationship  between  this  expansion  and  BWPT  and  also  introduce 
a non-iterative  method  to  obtain  Brillouin-Wigner  energies. 

In  chapter  4 we  propose  a method  closely  related  to  configuration  Interaction 
over  single  excitations  (Cl-S)  for  the  calculation  of  polarizabilities,  which  we  call 
Modified  CI-S.  This  techniqueallows  for  orbital  relaxation  through  the  inclusion  of 
single  excitations  on  a correlated  1st  order  Rayleigh-Schroedinger  wavefunciion. 


In  chapter  5 we  introduce  our  approximated  method  for  electronic  spec- 


iroscopy  that  we  call  Reduced-Expanded  Space.  We  try  to  incorporate  the  same 
level  of  correlation  captured  by  a ground  stale  1st  order  Rayleigh-Schrocdingcr 
wavefunction  to  the  single  excited  states. 

Finally,  in  chapter  6,  we  summarize  our  preliminary  results. 

In  1927,  Heitler  and  London  (24)  performed  the  first  quantum  mechanical 
calculation  on  a molecular  system,  namely  the  hydrogen  molecule.  This  paper 
marks  the  beginning  of  Quantum  Chemistry.  Theirs  was  a simple  valence-bond 
type  of  calculation. 

Much  of  what  is  done  in  this  thesis  though,  does  not  use  the  valence-bond 
approach  or  its  generalizations  [25,  26]  hut  is  rooted  in  the  Self-Consistent 
Molecular  Orbitals  method,  which  can  be  traced  to  the  ideas  of  Hartrec  and 
Fock  (27-30).  and  later  to  Roothaan  (31.  32).  and  the  molecular  orbital  picture 
of  Hund  (33)  and  Mulliken  (34-37). 

Actually,  the  first  quantum  mechanical  calculation  of  Heitler  and  London 
was  a correlated  calculation  as  defined  nowadays.  i.c.  correlation  energy  is  the 
difference  between  the  Hartree-Fock  limit  and  the  exact  nontelativistic  result.  But 
one  can  trace  to  Moller  and  Plesset  [38]  the  first  attempt  to  directly  obtain  the 
correlation  energy.  Perturbation  Theory  was  later  more  elegantly  developed  by 
Breuckner  and  Goldstone  [39.  40). 


Liiwdin’s  papers  [41.  42]  provided  an  elegani  way  lo  interpret  the  increas- 
ingly more  sophisticate  methods  like  Configuration  Interaction  (Cl)  [43-48],  Unre- 
stricted Hartree-Fock  (UHF)  (49).  Multi-Configuration  Self-Consistent-Field  (MC- 
SCF)  [50.  51],  Coupled-Cluster  Theory  (CC)  [52. 53]  and  methods  based  on  the 
Equations  Of  Motion  (EOM)  [54-57]. 

By  the  late  1960's  not  only  were  transition  energies  and  total  energies  being 
examined,  but  efforts  to  calculate  molecular  properties  had  intensified  [58-61]. 
Coupled-Perturbed  methods  based  on.  for  example.  Hartree-Fock  started  to  receive 
attention  [62-64], 

At  the  present,  the  effort  continues  to  obtain  more  accurate  properties  and 
Potential  Energy  Surfaces  for  larger  and  chemically  more  interesting  molecules. 

In  this  miroducuon,  I do  not  make  any  distinction  between  ab-initio  (all 
integrals  are  calculate  from  first  principles)  and  semi-empirical  calculations  (some, 
or  all  integrals  are  approximated),  because  in  this  thesis  1 consider  all  Hamiltonians 
in  second-  quantized  form  [65],  and  the  formalism  developed  are  identical. 

The  Electronic  Hamiltonian  from  the  Molecular  Hamiltonian 
The  Bom-Oppenheimer  Approximation 

The  Bom-Oppenheimer  Approximation  (BOA)  [66-68]  is  based  on  the  small 
mass  of  the  electron  relative  to  the  nuclear  mass,  which  allows  the  electrons 
to  adjust  instantaneously  to  the  nuclear  motion.  Hence  the  interaction  between 


^£(Ra  -Ra  );  A'  = l.A’-l  (5) 

is  defined  with  respect  to  the  center  of  mass  %k  of 


and  iik  >s  die  reduced  mass  of  nucleus  K and  the  total  mass  of  the  nuclei  already 


The  relative  coordinates  of  the  electrons  are  defined  with  respect  to  the  center 
of  mass  of  all  nuclei 

r,  = K,  - % (8) 


In  this  coordinate  system  the  Lagrangian  is 

-HW  m 

where  H and  r,  represents  collectively  all  nuclear  and  electronic  coordinates, 
respectively,  except  Rcm>  The  center  of  mass  motion  is  completely  separated 
since  it  does  not  appear  in  the  potential,  and  it  will  be  neglected  from  now  on. 
The  fourth  term  in  Eq.  9 is  called  the  mass  polarization  term. 


Introducing  the  conjugate 


We  construct  the  classical  Hamiltonian  by 


H = Y.PKitK  + Y.P^'-L  <12> 

explicitly 

+vyw  <°> 

where  M„„,  is  the  total  mass  of  all  nuclei. 

Replacing  the  momenta  by  their  quantum  mechanical  operators,  one  obtains 
the  quantum  mechanical  Hamiltonian  operator 


H = T,„  + «,,(/?)  + Hmp 


is  the  nuclear  Kinetic  Energy  operator; 

H„(R)  = -5  £>J  + V(ftr) 

is  the  Electronic  Hamiltonian  operator,  which  depends  parametrically  0 
nuclear  coordinates,  but  it  is  independent  of  the  nuclear  masses,  and 


(17) 


is  die  Mass  Polarization  operator. 

For  each  nuclear  geometry,  the  Electronic  Hamiltonian  possess  a complete  set 
of  eigenfunctions.  Hence,  one  can  write  the  solution  of  Schroedinger  equation 
H*  = E<b  (18) 

like  an  expansion  over  the  orthonormal  eigenfunctions  of  the  Electronic  Hamil- 

H,,(R)MR,r)  = WrWMIt,r)  OW 

where  the  expansion  coefficients  are  independent  of  the  electronic  coordinates. 
Different  basis  of  electronic  eigenfunctions  are  used  for  each  nuclear  configuration 
R.  One  should  be  careful  to  ensure  that  the  basis  functions  vary  smoothly  as  one 
changes  the  nuclear  configuration,  in  particular  when  degeneracies  occur  among 
the  eigenfunctions. 

In  principle  one  may  expand  the  total  molecular  wave  function  for  an  arbitrary 
nuclear  configuration  R in  the  complete  set  of  electronic  eigenfunctions  at  one 
particular  nuclear  configuration  R° 


*(*.r)=£CWlM*V) 


(21) 


When  the  summations  are  infinite,  both  expansions  (20  and  21)  are  exact.  If 

flexible  over  the  complete  PES,  because  of  the  extra  freedom  of  having  a different 
electronic  wave  function  for  each  nuclear  configuration. 

To  determine  the  expansion  coefficients,  following  Bom.  we  insert  the  expan- 
sion Eq.  (20)  into  the  Schroedinger  equation.  Eq.  (18).  multiply  on  the  left  by 
r)  and  integrate  over  the  electronic  coordinates.  The  resulting  differential 
equations  for  the  expansion  coefficients  arc 

[- 257  E v*  + WW  + j f7cw<*)  - £]  '*(«)  = ~jj  E <VWx,(«) 

' (22) 

C„.(R)  = Y.  [-'^(«)-VA-  + + D„.(R)  (23) 

coupling  matrix  elements  (NACME's) 

= (ti„(fl.r)|  - iVA-Kv(/?,r))  (24) 

the  second-order  NACME's 


(25) 


polarization 


°£™ --jE-t'MA'II 


The  above  form  of  the  Schroedinger  equation  is  exact.  No  approximations 
have  yet  been  introduced. 

The  Schroedinger  equation  is  greatly  simplified  if  we  neglect  all  matrix 
elements  This  results  in  the  Bom-Oppenheimcr  adiabatic  wave  function 

r)  (27) 

where  the  nuclear  wave  functions  v/»(^)  “re  the  solutions  of  the  equation 

H?°WxiAR)  = (28) 

= -^jE  V*.  + Wn(fl)  (29) 

For  each  electronic  state  there  is  a set  of  nuclear  wave  functions  representing 
different  nuclear  states.  These  are  enumerated  by  the  index  /. 

In  the  Bom-Oppenheimer  approximation  the  electronic  energy  acts 

as  a potential  for  the  nuclei.  This  is  equivalent  to  assuming  that  the  electrons 
adjust  instantaneously  to  the  motion  of  the  nuclei.  There  is  thus  a different 


Bom-Oppenheimer  Hamiltonian  for  each  electronic  state.  The  diagonal  elements 
C pul  ft  ) are  sometimes  retained  in  the  Bom-Oppenheimer  Hamiltonian  since  they 
do  not  couple  different  electronic  states,  and  this  is  referred  to  as  the  Adiabatic 
Approximation.  Note  that  in  the  Adiabatic  Approximation,  the  potential  for 
the  nuclei  depends  on  the  nuclear  masses  since  the  correction  C„r(R)  is  mass 
dependent. 

Second  Quantization 

Second  Quantization  uses  the  occupation-number  representation.  It  is  a very 
convenient  representation  to  present  most  of  Quantum  Chemistry,  and  provides 
an  unifying  formulation  for  both  semi-empirical  and  ab-initio  methods. 

In  this  section,  we  briefly  introduce  the  essential  ideas  of  second  quantization 
as  well  as  define  our  notation.  A more  complete  introduction  to  second  quantized 
methods  can  be  found  elsewhere  [57,  65.  68.  69J 

In  first  quantization,  the  electronic  Hamiltonian  is  written  as 

n.,w = E *«) + Y.9VJ)  + E (30) 

where  the  last  term  is  just  a constant  for  each  nuclear  configuration  and  can  be 
added  at  the  end  of  the  calculation.  Hence  we  are  interested  in 

(31) 


h = E mo + 


on 

is  an  one-electron  operators  and 

<33> 

is  the  two-electron  part  of  the  Hamiltonian. 

Electrons  are  fermions,  so  a wave  function  describing  their  motion  must  be 
antisymmetric  under  interchange  of  any  two  panicles.  This  is  a property  of 
functions  constructed  as  a linear  combination  of  Slater  determinants 

x.HEWn.x, x„)  (34) 

where  x,  is  the  coordinate  (spacial  and  spin)  of  the  i-lh  electron  and 

tMx  1.*I *■)  = (A'!)-1/!dct  |<S|(Xl)  . . .<S,(x„-l)$r(x„)|  = Idi . - -0,6,\ 

(35) 

The  {o.}  is  an  orthonormal  (without  loss  of  generality)  one  electron  basis  set. 
These  one-electron  function  are  called  spin-orbitals. 

Slater  delerminantal  wave  functions  involving  orthonormal  spin-orbitals  o, 
can  be  represented  in  terms  of  products  of  creation  operators  on  the  so-called 
vacuum  ket  |uac). 

rV...f*|uoe)  <-*  |*  . - 0,0, | 


(36) 


The  fundamentals  of  the  second  quantization  are  summarized  in  table  1. 


■nible  1: 

mndamentals  of  second  quantization 

Basis  vectors 

|n)  = |ni,n2,...,nm);  n,  = 0,1  for 
i — 1,2 m 

Inner  product 

(n|k)  = t..h 

Creation  operators 

nm)  = 

r(n)|n, 1 nm) 

T(n)  = (—1  1 

•'|n ,li,...,nm)  =0 

Annihilation  operators 

HnJIn, 0 nm) 

i|nt 0 nm)  = 0 

Anticommutation  relation 

= ‘'j  + J>’  = &ij 

= i# + J>  = 0 

Number  operators 

i»i|n)  = n,|n) 
Ei'i|n)  = E".ln) 

Vacuum  state 

(!>tic|uoc)  = 5 

The  second-  quantized  version  of  any  operator  is  obtained  by  simply  de- 
manding that  the  operator,  when  •'sandwiched"  or  "bra-ket-ed''  between  bra  and 
ket  vectors  of  the  form  J] (■>"!>'  = (II  r'li’ac))'  and  I]  r'ltiac),  yields  exactly  the 


same  result  as  arises  in  using  (he  first  quantized  operator  between  corresponding 
Slater  determinant  wave  function.  The  results  are  shown  in  Table  2. 


Table  2:  Characteristics  for  operators  in  first  and  second  quantization 


First  Quantization 

Second  Quantization 

one-  electron  operator 

one-  electron  operator 

£/(*.*)  = £/(i) 

E(*r.l/l*,)pV=E/„P,9 

two-  electron  operator 

two-  electron  operator 

iEs(r,,»i;rJ,e,J  = {EjW) 

i E<(®p<i«l9l«W<)pVsr  = 

JpE  Snr.p'q'sr 

operators  are  independent  of  spin- 

operators  depend  on  spin-  orbitals 

orbitals  basis 

basis 

operators  depend  on  the  number  of 

operators  arc  independent  of  the 

electrons 

number  of  electrons 

exact  operators 

projected  operators 

So  the  second-quantized  form  of  the  Electronic  Hamiltonian  (except  for  the 
constant  nuclear  repulsion  term)  is 


H = h„p'q  + -(w||rs)pVs 


(37) 


*)*<*■> 


dx,  (38) 


(wlM  = / «;(X1  W(x,)j^(l  - «,)«(*, )«(*,)  *.*C2 


CHAPTER  2 

METHODS  OF  APPROXIMATION 


The  simple  pictures  used  by  chemists  and  physicists  are  based  in  the 
independent-particle  model.  The  most  widely  used  independent-particle  model 
is  the  Hartree-Fock  method.  It  generates  the  "best"  spin-orbitals  set. 

The  Hartree-Fock  method  is  important  not  only  because  it  is  in  itself  a good 
approximation  but  also  because  it  is  the  starting  point  of  many  more  accurate 

Here  we  will  briefly  describe  the  Haruee  Fock  method,  and  review  some 
of  the  most  successful  past  Harliec-Fock  methods  for  including  the  correlated 
motions  of  electrons.  1 focus  on  the  Configuration  Interaction  Method  (matrix 
diagonahzauon),  Perturbation  Theory,  and  the  Coupled-Cluster  Ansatz. 

symmetry  Adapted  Coupled-cluster — Configuration  Interaction  (SAC-CI)  [1-18]. 
where  we  proposed  some  corrections  to  the  original  work  of  Nakalsuji  and  co- 
workers. We  also  introduce  a diagrammatic  way  to  obtain  the  SAC-CI  ground 
state  direedy  over  molecular  orbitals. 

Hartree-Fock 


[27-32,  65,  68-71] 
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of  spin-orbitals  for  which  the  energy  is  an  extreme,  when  the  wavefunedon  is 
approximated  by  a single  Slater  determinant  constructed  with  those  spin-orbitals, 
subject  to  the  orthonormality  condidon. 

One  writes  a trial  wavci  uncuon  as 

$ = Iti'h  (M4j)  = fij  (41) 

which  assures  that  = 1.  By  applying  the  Eckart's  Variational  Principle  (VP) 
[72-78]  with  the  reslricdon  that  the  spin-orbitals  conunue  to  be  an  onhonormal 
set.  i.e.  (o,|Cj)  = (t\j)  = 6tJ,  one  obutins  the  Hartree-Fock  equations.  This  is 
done  by  constructing  a Lagrangian  and  introducing  Lagrange's  muldpliers  { ear } . 
explicitly 

= -][>,(* |/>  (42) 

H{v]  = (ti|  H |t >)  = »**  + j(«| \kl)  (43) 

6C  = 6H  (44) 

If  one  constrains  6H  to  be  real,  then  i‘w  = in.  So  varying  the  orbitals  and 
setung  6C  = 6H  = 0,  one  geus 
0 = (<p|{  [*  + (/ rfx,  *;(x2)^(l  - f*,2)d,(*a))]  Ip)  - 

(45) 
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Since  the  variations  (££|  are  unrestricted  then 

/(xi)d,(x,)  = £„d,(x,)  (46) 

where  f{x)  is  called  the  Fock  operator  and  is  defined  as 

/(xi)«M*i)  = [h(xi)+  (J dxi  «s;(*2)^(l  -/’12)d,(x2))]dp(x1)  (47) 

For  the  closed-shell  restricted  Haruee-Fock  (RHF)  case  or  for  the  unrestricted 
Hartree-Fock  (UHF)  case  i is  an  hermitian  matrix,  it  can  be  diagonalized  by  a 
unitary  transformation  (e  = UsU1 ; UU1  = U'U  = 1).  resulting  in  the  canonical 
Haruee-Fock  equations 

/(*i)Oj(Ji)  = £/  Pl{T,)iu  (48) 

A similar  equation  can  also  be  derived  for  the  general  open-shell  restricted 
Hartree-Fock  case,  but  this  derivation  is  more  complex  as  is  the  form  of  /(x] ). 
For  systems  with  an  even  number  of  electrons  (A'),  the  RHF  assumes 

|0>  = tor  = |df  dfd . . . dSyX/jl  (49) 

and  the  energy  is  then  given  by 


EHF  = (0|tf|0)  = eu  - ±(M||M) 


ZkZl 

Rkl 


(50) 


«M-*»  + (H||H)  (51) 

The  exact  solution  of  the  RHF  equations,  i,e.  o and  B spin-orbitals  have  the 
same  spacial  part  (orbital),  result  in  the  Hartrcc-Fock  limit  energy.  The  difference 
between  the  RHF  limit  value  and  the  exact  value  of  the  nonrelativistic  Electronic 
Hamiltonian  is  called  correlation  energy  [42],  So,  by  definition.  Restricted 
Hartree-Fock  energy  has  no  correlation,  although,  physically  it  correlates  the 
motion  of  electrons  with  different  spins. 

The  solutions  of  the  Hartree-Fock  equations  are  usually  attained  by  an  expan- 
sion of  the  molecular  orbitals  (o,)  over  atomic  orbitals  explicitly 

0 = x-C  ; = X,C,P  (52) 

With  this  approximation  the  Fock  equation  for  the  closed  shell  case  is  written 

FC  = eSC  (53) 

F„  = (Xpl/IXs)  = <XplMX,)  + £ C;„C„<x,X.||XA>>  (54) 


S„  = (Xpl.Xfl) 


(55) 
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If  all  the  integrals  of  equation  54  are  evaluated  from  first  principles,  it 
is  an  ab-initio  method  [79).  Approximations  are  made  to  increase  the  speed 
and  accuracy  of  the  method  and  such  approximations  to  the  one-electron  part, 
Dra  = YL  C‘VC,V  and  the  two-electron  integrals  results  in  a scries  of  semi- 
empirical  methods  [80-83]. 


Configuration  Interaction 

For  most  situations  the  HF  method  is  able  to  give  qualitative  agreement  with 
experiment.  Examples  are  geometric  structure,  vibrational  spectroscopy  and  the 
energy  differences  connected  with  conformational  changes  such  as  inversion  and 
rotational  barriers.  Proton  affinities  and  hydration  energies  are  also  well  described 
by  the  Hartrec-Fock  method.  Common  to  all  cases  for  which  the  HF  method  is 
applicable  is  that  no  drastic  rearrangement  takes  place,  while  the  system  changes 
from  state  A to  state  B.  No  electron  pairs  arc  broken  or  formed,  and  the  molecular 
orbitals  are  only  slightly  modified.  The  error  connected  with  the  HF  approximation 
remains  almost  constant,  and  a qualitative  agreement  can  be  achieved.  Many 
interesting  problems,  however,  do  not  belong  to  this  category,  e.g.  ionization 
energies,  electronic  excitation  energies  and  especially  important  chemical  reactions 
involving  the  formation  of  covalent  bonds. 


arc  then  necessary.  Historically,  Configuration  Inter- 


action (CI)[41,  43,  45,  65.  68-71]  was  one  of  the  first  correlated  method  based 
in  the  molecular  orbital  theory,  and,  in  fact,  was  originally  the  title  given  to  all 
post  Hartree-Fock  methods.  Today,  Cl  usually  refers  to  methods  of  the  variational 
type  that  results  in  matrix  diagonalization. 

Expanding  the  wavefunclion  as 

Id.)  = (Co  + C?o'r  + Ctfa'b'ji  + . . ,)|0)  (56) 


where  (C.)  are  variational  parameters,  and  applying  the  variational  principle  to 


(57) 


yields  the  secular  equation 


where  t’  is  the  row  vectors  of  configurations  and  C is  the  column  vector  of 
variational  parameters.  The  Cl  wavcfunction  can  be  written  as 

* = vC  (59) 

while  H is  the  Hamiltonian  matrix  between  the  different  configurations 


H„  = Iv*,} 


(60) 
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The  Cl  wavefunction  is  size  extensive,  only  when  all  configurations  within  a 
basis  set  are  included,  which  is  called  Full  Cl  (CIF).  The  eigenvalues  of  the  secular 
equation  are  upper  bounds  of  the  electronic  states  order  by  order  [43, 44,  48). 
Size  cxlensivity  refers  to  the  property  that  A'  non-interacting  subsystems 
HaW 'a)  = Ea\<I>a) 


(61) 


HkI’I’k)  = En\<Pn) 

should  have  total  energy  given  by  the  sum  of  its  pans,  i.e„ 


(Ha  + Hb  + ...  + Hum  = {Ea  + £fl  + . . . + (62) 


Perturbation  Theory 

Perturbation  Theory  (PT)  [38-40,  47.  65.  68.  69.  84-89]  assumes  that  the 
initial  guess  for  the  wavefunction  is  indeed  a good  one.  and  expands  the  energy 
in  terms  of  correction  to  this  initial  guess. 

Rewriting  the  Scrocdinger  equation  as 

(C-H,M  = (V-E  + eM)  (63) 

where  f is  a shift  parameter  and  assuming  intermediate  normalization  , i.e. 


W«|0)  + W> 


(64) 


<3|0)  = 0,  and  (0|i b)  = (0|0)  - 1;  one  obtains: 


<?|ti-)  = (£-ff.)-’Q(V-£  + £)W  (65) 

A£  = £-£„  = (Oil'll)  (66) 

£»  = (0\H\0)  (67) 

Q = \-  P = l- 10)(0|  (68) 

(e  - Ho)~'Q  = R is  defined  as  ihc  Resolvent  . yielding 

W>-  fl(?(V-£  + £)|tir>  (69) 

For  Rayleigh-Schroedinger  Perturbation  Theory  (RS)  e = £0.  thus 

Rss  = {Eo-H,r1Q  (70) 

For  BriUouin-Wigner  Perturbafion  Theory  (BW)  c = E,  thus 

Raw  = flatv(£)  = (£- tf.r1  <2  (71) 


Iterating  over  the  wavefunction  equation  69  yields 


W = |0)  + flt?(V-£  + e)|0)  + ...  (72) 

= plRQ(V-E+e)m  (73) 

= f>,*,>  (74> 

AE  = (OU'lf)  (75) 

= f>|V'  |fl(?(V-  £ + £)]"  |0)  (76) 

= f 'J"  (77) 


Since  the  HF  wavefunction  in  many  cases  provides  a reasonable  description 
of  the  electronic  structure,  it  is  natural  to  improve  on  it  by  perturbation  theory. 
In  Mpller-Plesset  and  Many-Body  Perturbation  Theory  (MBPT).  the  partitioning 
of  the  Hamiltonian  lakes  this  in  consideration,  explicidy 


h!"’  = /«/>'<) 


(78) 
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VUP  = j(m||m)p',V  - {pk\\qk)p'q  (79) 

E'mp  = -\(klW)  (SO) 

Ehf  = E°mp  + E'mp  (81) 


Rayleigh-Schroedingcr  perturbation  theory  can  be  shown  to  be  size-extensive 
order  by  order  [68).  but  it  is  not  an  upper-bound.  Brillouin-Wigner  Perturbation 
Theory  is  not  size-extensive  but  its  odd-orders  are  upper  bound  for  the  energy. 
For  example,  equation  81  implies  that  Hartree-Fock  is  size-extensive. 

The  main  disadvantage  of  perturbation  theory  is  the  assumption  of  a Hamil- 
tonian partitioning.  + V.  that  is  rapidly  convergent.  To  insure  this  may  require 
a reference  function  of  several  configurations  and  such  theories,  at  higher  orders, 
become  quite  complex. 

Coupled  Cluster  Theory 

The  basic  idea  of  Coupled  Cluster  theory  [52.  53.  65.  68.  69.  89,  90)  is  to 
treat  the  interaction  among  the  electrons  as  an  interaction  among  electrons  within 


clusters  and  a coupling  between  these  clusters.  The  coupling  between  the  clusters 
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of  electrons  is  obtained  by  writing  the  wavefunction  as  an  exponential  of  the 
cluster  operator  T acting  on  a single  determinant  reference  wave  function.  The 
reference  wave  function  is  most  often  assumed  to  be  a Hartree-Fock  state 

\CC)  = exP(T)|0)  (82) 

where  the  cluster  operator  T consists  of  one-electron,  two-electron A'-electron 

clusters 

T = T,+T2  + ...  + Tk  (S3) 

T,  = t’a'i  (84) 

T.=  '-  aWj  (85) 

Again  the  labels  i,j.  Is,  I and  o,  t,  c,  il  are  used  for  occupied  (hole)  and  virtual 
(particle)  spin-orbitals  in  the  reference  determinant,  respectively.  The  cluster 
operator  truncates  at  excitation  level  A',  the  number  of  electrons  in  the  system. 
In  a shorthand  notation  the  cluster  operator  can  be  written  as 

r = ^r,  = y^f,T,  = tr  (86) 

where  (j  denotes  the  cluster  amplitudes  at  excitation  level  i and  t.  is  the  corre- 
sponding excitation  operator.  The  excitation  operators  satisfy 


(87) 


(O|r*T„|0)  = t. 


(88) 


|r(l,T„)  = 0 =■  [r„7'„]  = 0 (89) 

Hence  ihe  Schroedinger  equation  takes  the  form 
«e*p(7-)|0)  = £exp(7-)|0) 
or  rewriting  in  a more  convenient  form 

exp(-r)//exp(7-)|0)  = £|0) 

which  defines  an  effective  Hamiltonian 

H,„  = exp{—T)H  exp(T)  = « + [//,7')  + i[[//,r],r)+ 

i|p,r].ri,ri  + j;||[[//,ri,r],r],r] 

The  truncation  after  four  commutators  is  exact  regardless  of  the  excitation 
level  at  which  T is  truncate.  The  truncation  is  due  to  the  fact  that  the  Hamiltonian 
contains  at  most  a two-electron  operator. 

The  Coupled  Cluster  equation  is  solved,  in  general  by  requiring  that  Eq.  91 
holds  in  the  K-fold  excitation  manifold,  where  K is  the  excitation  level.  Explicitly 


(90) 


(91) 


(92) 


E = <0|  exp  (-T)H  exp  (T)|0>  = (0\H  exp  (T)|0) 


(93) 


(/i  | exp  [—T)H  exp  (T)|0)  = 0 


(94) 


W = (0|r* 


(95) 


Hie  energy  in  the  Coupled  Cluster  model  is  nol  variational  in  the  sense  that  it 
is  not  an  expectation  value  of  the  Hamiltonian.  The  minimization  of  the  quantity 
£ _ (0|  exp  (r)*//  exp  (r)|0)  (96) 

(0|  exp  (T)*  exp  (T)|0) 

would  give  rise  to  a variational  energy  with  the  CC  ansalz  of  the  wavefunclion. 

The  CC  equations  are  non-linear,  hence  they  have  more  than  one  solution. 
There  is  at  least  one  solution  which  is  size-extensive,  which  is  the  solution 
obtained  by  the  usual  techniques.  This  solution  can  be  viewed  as  an  infinite 
summation  of  MBPT  series  [68]. 

The  advantage  of  the  truncated  CC  approach  over  the  truncated  Cl  method  is 
that  the  former  is  size-extensive  and  give  a very  compact  parametrization.  How- 
ever. there  is  not  a complete  satisfactory  formulation  of  the  CC  approach  capable 
of  treating  excited  slates  or  near-degeneracy  situations  [68].  The  advantage  of 
the  truncated  Cl  method  is  that  it  always  supplies  a well  behaved  upper  bound 
to  the  energy. 
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Spin-Symmetry  Adapted  Coupled  Cluster 

One  of  the  more  promising  CC  based  methods  for  excited  states  is  the 
Spin-symmetry  Adapted  Coupled-cluster  - Configuration  Interaction  (SAC-CI) 
approach  [1-18].  We  will  concentrate  in  the  singlet  case  here,  and  present  a 
diagrammatic  representation  of  the  SAC-CI  method. 

The  SAC-CI  method  is  used  to  calculate  transition  energies  which  takes 
advantage  of  the  spin  symmetry  to  reduce  the  number  of  integrals. 

The  SAC-CI  method  was  introduced  by  H.  Nakatsuji  [1-18].  It  has  a Coupled 
Cluster  wavefunction  for  the  ground  stale.  The  muliielcctron  basis  set  is  generated 
from  single  excitations.  The  excited  states  are  obtained  in  a fashion  paralleling 
the  Cl  method  [7], 

The  ground  state  wavefuncuon  is  written  as 

|v'-  ) = 6sexp(S)|0)  = 6il/)  (97) 

where  S = St  4 Si  4 . . . + Sjv  is  the  spin-symmetry  adapted  cluster  operator  with 
Si  = CfS‘K  Si  = CflS"}',  etc;  6,  is  the  spin-projector  of  the  ground  state; 
and  |0)  is  the  reference  single  determinant. 

S“'  is  S"1  = 4 «Jig)  for  singlet  excitations.  S°'  = a'jg  for  triplet 

excitations,  etc. 


(H-E)0  exp(S)|0)  = 0 (98) 


One  solves 


30 

10  obiain  the  spin  adapted  CC  amplitudes  (C)  by  projecting  the  above  expression 
on  a set  of  lower  order  excitations  out  of  the  reference  function,  as  in  equations 
93  and  94. 

The  multieleciron  basis  set  {d,}  for  the  excited  (or  ionic)  states  is  generated 
by 

|di)  = AtQ6rfl*|d-')  (99) 

where  A t is  a normalizing  constant,  Q (Q  = 1 — K;9)(t^l)  projects  into  the 
orthogonal  space  to  |t!>  ),  d,  is  the  projector  into  the  desired  spin  symmetry, 
and  R J is  the  symmetry-adapted  excitation  (or  ionization  or  electron  attachment) 
operator. 

Thus,  one  can  write  a wavefunction  for  an  excited  state  as 

\<i’c)  = 1 9k)dk  + ity</e  OOO) 

d =0  if  the  excited  state  and  the  ground  state  do  not  have  the  same  spin 
symmetry  or  the  ground  slate  was  obtained  variationally. 

In  any  case,  one  has  to  diagonalize 

Hd  = d£  (101) 

where  H = {o,\H\o  ).  and  d = (dt,d.,.  - ,dm).  to  obtain  the  excited  state 
energies  and  wave  functions.  Such  a procedure  re-inuoduccs  unlinked  clusters, 
and  the  method  looses  its  size-extensivity,  although  this  appears  to  be  of  minor 


importance  in  numerical  examples.  However,  diagonalization  appears  to  be  a 
necessary  step  to  obtain  non-interacting  spectroscopic  states. 

For  a singlet  ground  state  and  singlet  excited  stales.  Nakatsuji  [1-4]  assumes 
the  following 

) = exp(S)|0}«(l  + $ + 5[S2)|0)  002) 

where  5 at  S\  + 52- 

5,  = C?S"\  and  S2  = and 


(w.(.+er«,+(ef+jacM,sr+  (1M) 

CfCj|5r,S*,Sjt  + ^Cjfsr's*)srs'l)|°> 


(105) 


To  obtain  the  cluster  amplitudes,  he  further  requires  tha 
(0|W  - E|d>,)  = 0 
(Sf'0|H  - Bid-,)  = 0 
(Sj’sJ*0|A  - B|d-,)  = 0 

which  means  requiring  that  (H  - E)|d-,)  = 0 in  the  sub-space  spanned  by  the 
reference  function,  single  and  double  excitations  with  respect  to  this  function. 


The  amplitudes  can  also  be  obtained  variationally,  by  applying  the  variational 


principle  to 


W«C})  = 


(106) 


The  singlet  excited  states  basis  set  is  generated  by 


|/)  = A.o<?S°Vj) 


(107) 


Hence.  |ti'e)  = |o“  )d“  + K'.)d0.  The  solutions  are  obtained  by  diagonalizing 
Hd  = d£.  If  |«i  ) is  obtained  variationally.  then  </„  is  considered  zero  because 
it  will  not  couple  with  the  other  functions  through  the  Hamiltonian. 

The  Nakatsuji  ansatz.  for  singlet  ground  stale  |t’s)  = exp( S)|0)  ss  (1  + S + 
4iS2)|0)  leads  to  disconnected  terms  in  the  reference  energy  expression,  since  the 
exponential  operator  is  truncated  at  this  order. 

1 suggest  that 


|,y  = expose  (l  + S + ls2  + 


(108) 


ltil)  = {i  + c“s“,  + (c“‘  + 5c“c‘)s'1s‘,+ 


(109) 
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if  we  assume  S re  S,  + S2.  This  new  expression  for  ihe  ground  slate  leads  to 
connected  terms  only,  in  the  coupled  equations.  That  is,  all  the  unlinked  terms 


cancel  in  lower  order. 

(0|tf  - £|0f)  = 0 (110) 

(S°'0|/y  - £10,)  = 0 (111) 

(S1°,Sj,0|/f  - £|v''j)  = 0 (112) 

and  the  resultant  energy  is  size-extensive,  at  least  until  matrix  diagonalization. 

In  terms  of  the  molecular  orbitals  we  have 

f„  = h„  + ((pi|f>))  <ll3) 

((TC|rs))  = 2(pg||rs)  - (pg||sr)  (114) 

E"  = Ehf  = (0|H|0)  = 2/„  - ((V|ij))  (H5) 

(pr|gs)  = Ws)  = J *!(iradj(r,)dj(r2)-|-d,(ri)<i,(r2)  (116) 


Assuming  Eq.  109  and  using  Eqs.  110.111,112.  one  obiains  a sol  of  coupled 
equations  in  terms  of  molecular  orbitals  instead  of  spin-orbitals 

1. {0|H  - E|v'’g>  = 0 

£ = E°  + V5C°/ . + (Cat  + 5C“c‘)((ii|u4)) 

1 V i » 1 (1,7) 

E = 2/..  - ((ii|y))  + V5C*/ja  + (C ..  + -C  .C.){{ij\at>)) 

2. <0|S'(H-E)|t!'g)  = 0 

s/2/  , + Caf,  -Cff.  + C°((i/|om))+ 
mf  m fa  i im  i 


+C^(0|H-£|i'j)=0 


(118) 


3.{0|S^S‘(H-E)|l'g)  = 0 


((/9|m")l+ 

+^{C°((/9|«n)>  + C“(</s|m0,))+ 
--L(C{  ((.'slmn))  + C9.  ((i/|nm)))+ 


(119) 


+{-2C*  + C®7.}/  ,n  + {-1C!.an  + C°[}fim+ 

-{"2Cmn  +Cml)/  o “ <_2Cmn  + Cmn,//a+ 

+{2c'“  -C^-  5C  “ }((ij|«n))+ 

+{2C®“  -C“®  -ic“c®}((.7|am))+ 
-C£(W|«))  - {C°®.  + 5C  “ C®}((i/|n0))+ 
-c'“((ij|om))  - [C‘nf.  +\cBnCf.  )((ij|m«))+ 

+{c“‘n  + ic“c‘}((/s|0i))+ 

+{C,a.  +\cf.C3.m\mn))+ 

+4={— 2 CaC9f  + CaCf9.)f.  + 

71  mm  m m ifl 

+-^{-2C“c^®.  + C“c^.}/.o  + 

+-^{-2C®C^°  +C^C®°  )/.o+ 

+-^{-2C^C°®  +CSiC„n)fia+ 


^{-2Cf.CS“.  +C®C^“}((.j|«m))+ 
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+^g{-at7*c*  + c*c^K(«|ow))+ 

+ +c“c3/.wi\am))+ 

+-^{-2  C'.C^.  4 C^C*a.}((ij|an))+ 

+4={-2  CaC9f  +C°C^}((i;|an))+  (124) 

y/2'  > jm  I jm 

+-j={C9.CaJ.  4 CfiCj?n+C^Cf9)[{ij\an))+ 

--^<-2C  ‘ C“®  + C *cf®  )((i/|«&))+ 
--L{-2C“Ci^|+C“Cn^i)((.7l«6))+  (125) 

__>_{C*C  ^ + C ‘ c”*  4 C“c^.}((i/|at))+ 

--^{-2 C*c”£  4 C * c”^}((ij|«i-))4 
--^{-2C“C^  4C"c^/n)((isl<i4))4 
--^(c'c  4 <?“  4 CbnCaJ.mig\cb))+ 

+J^{cV*c(((y|«m))4C“c{c*((ii|«n)))4 
_2^51CiCmCn(('/|0<,)>  + C iCmCn((,9l<,i)))+ 


+ 1_SCm71Ci{  +L'™«1'  ij  ' 

-<-Ki'cS+oScSll®<‘»+ 


+ii-!c'c'c  » + c-c>c»  +c‘yA‘  !«‘W+ 
+lt— ic-c V'  + c'c*c  *'  4 c'.c’.c^ KWMI+ 
+i(-w  vV*  4cv»cj* +c;vV'K(«i.t»4 

+I(-w‘o*o^4i:io'oJ»4cio>c;'](((jW)4 
+Il  4 CVjc  ’*  4 OV’c“  )«ljl-‘))4 

+I{-K;c"C*40;C*C^4cVfc"K(i/l*»))4 
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+jC|c9cV‘((0'|<i6))+ 

+2 {C^  + 5^C®}(0|tf  - E\^g)+ 

+C^<0|S*(«-£)W>fl)+ 

+C9n{0\sjn{H-E)\>!’  )+ 

-!<(<(" -£>|^ )+ 

" 5Cm(0|Sn(A"£)lV' 

Diagrammatic  Rule  for  a SAC -Cl  Singlet  Ground  State 

The  expressions  in  Eqs.  117.118,119  can  be  derived  using  diagrammatic 
rules,  that  we  introduce  here.  The  rules  are  different  from  the  usual  ones  in 
PT,  because  we  are  using  molecular  orbitals  instead  of  spin-orbitals,  and  have 
spin-eigenstates  restrictions. 

Diagrammatic  Rules: 

1.  Full  arrows  indicates  sum  over  index 

2.  Sign:  (-iy*Ha*Po  where 


is  the  number  of  holes  summed 


p is  the  number  of 
open  lines) 


lines  divided  by  2 (or  the  number  of  pairs  of 


1q  is  the  number  of  independent  a loops 
Pa  is  the  number  of  independent  a paths 

3,  and  open  lines  are  never  equivalent 

4.  (2)Iq  *po  ; all  continuous  paths  or  loops  can  be  o or  0,  except  those  where 
the  panicle  and  respective  hole  are  not  in  same  path  (or  loop) 

5-  (^)p  : each  pair  of  open  lines  contribute  with  this  because  S/, 1 = 

£(/!». + />*»#> 

6.  T contributes  with  (tn|/'|ouf)'  = (ou/|F|m) 

7.  j conuibutes  with  |((ouf  |tn))  = £((in|ouf))' 

8.  C,  contributes  with  while  C , contributes  with  1 


In  the  following  examples,  1 will  clarify,  hopefully,  the  meaning  of  these  rules. 
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Figure  1:  Diagrams  over  Molecular  Orbitals:  example  1. 

1-  Sign:  (-l)h*P*|a*',Q=  (-l)l*0',*0=  + 

2-  No  equivalent  lines  = 1 

3-  (2)'Q-PQ  = (2)'*°=  2 

4-  No  open  lines  = 1 


Hence  this  diagram  represents  2 /,, 


Figure  2:  Diagrams  over  Molecular  Orbilals:  eicample  2. 


1-  Sign:  (-l)"*P*'o*PQ=  (-1)2*°*J*°=  + 


3.  (2)'q-Po=(2)2'0 


6-  "j 


Hence  Ihis  diagram  represents 

+(-^)2  x (2)1  x 1 x i((ij'|«4))CJ?  = ((y|«»))Cf/  033) 
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Figure  3:  Diagrams  over  Molecular  Orbiials:  example  3. 


The  above  diagrams  clarify  Ihe  concept  of  independent  paths.  One  has  to 
assign  connected  pairs  such  as  (m.  0.  (n.g)  for  all  double  excitations  associated 
with  <0|S's;(ff  - £)|ii>j)  = 0.  For  the  left  diagram,  the  (m,f)  and  (n.g)  pairs 
can  be  a or  B spin  independendy,  two  independent  o paths.  The  right  diagram 
has  only  one  independent  o path,  because  if  (m.O  is  o spin,  then  so  is  (n.g). 

So  the  left  diagram  yields 


(-1)' 


1+2+0+2  1 0.  2 + 0 
I 1^1  (4) 


(134) 


and  the  right  diagram  yields 


(-1) 


CHAPTER  3 

A PERTURBATION- VARIATIONAL  METHOD  FOR  THE  GROUND  STATE 


Introduction 

The  ground  state  potential  energy  surface  is  and  has  been  subject  of  intense 
research  over  the  years.  A reliable  description  of  the  ground  state  at  all  inter- 
nuclear  distances  is  important  to  describe  equilibrium  geometries,  rovibrational 
spectra,  and  the  details  of  reactions. 

In  this  chapter  we  revisit  perturbation-variational  theories  (19.  91-94.  20. 
22.  23).  and  Brillouin-Wigner  perturbation  theory  [84-86).  that  are  capable  of 
yielding  a large  fraction  of  the  correlation  energy  rapidly  and  check  the  reliability 
of  these  methods  in  yielding  the  ground  state  correlation  energy.  Our  extended 
goal  in  examining  these  methods  is  to  obtain  a balanced  theory  for  molecular 
spectroscopy,  and  this  influences  some  of  the  ideas  we  present  here. 

In  section  2,  we  present  a non-iterative  procedure  for  Brillouin-Wigner  pertur- 
bation theory.  In  section  3.  we  analyze  a variational  method  using  a Krylov-type 
expansion  and  its  connection  to  perturbation-variational  procedures.  In  section  4. 
we  present  some  of  our  results.  Finally,  in  section  5.  we  summarize  our  results 


and  present  our  conclusions. 


Brillouin-Wigner  Penurbalion  Theory 


Remembering  that  the  Resolvent  is  defined  as 

(e  = fl  <I36> 


For  Rayleigh-Schroedinger  Perturbation  Theory  (RS)  £ - Ec.  and 


R„s  = (E.-H.r'Q  <137> 


For  Brillouin-Wigner  Perturbation  Theory  (BW)  e = E,  and 


Rbw  = Rbw(E)  = IE-H,)-‘Q 


(138) 


Remembering  that 

M)  = = f>  Q (V  ■ - E + £)]"|0)  (139) 


AE  = f)£"  = (0|VW 

»=i  (140) 

= f>|V  [RQ(V-E  + e))’\0) 
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Assuming  lhat  A E is  small  by  carefully  partitioning  the  Hamiltonian,  which 


= E(-A  £)"+,(«fls)"+J 


Raw  = 53(-A  £)"-’(flRs)"  (142) 

law"')  = (flflH.V')W|0)  = ^f;(-A£r‘(fl*i)"  j v|  |0)  (143) 

and 

£fliv  = (0|V|BH'w)  = (0|V(JfcwV)»|0)  = (0|V  m-ABr-'^rV  |0) 


(144) 


E'bw  = E'rs  = (0|V|0)  (145) 

E\w  = fjt-D-lfiV  + E%w)’-' '(OlVffiflsJ-VIO)  (146) 

Ebw  = E E(-l)"+m-J(£V+£!tv+4w),+m_J(0|l'(fl«s)nV(iiflsri'|0) 


(147) 


e*bw  = f;  {-*E)l+‘+n,-3(Q\V(RKs)tV(RKs)‘V(RnsrV\0) 

t.i.l-l.m-1  (148) 


Ej,s  = (0|1'««SV|0)  (149) 


E3„s  = (0|Vflj,sV7?fls('|0)  - £Is(0|K(*im),K|0)  (ISO) 


EtRS  = {0\VI{Rs{V-{V))RRs(V-{V))RKSV\0)-{0\VRRsV\0mV(ltKs)1VHl) 

USD 


From  equations  143  and  144,  one  can  see  that  the  evaluation  of  the  Brillouin- 
Wigner  terms  that  appear  from  the  use  of  expansion  142  are  as  difficult  as 
the  corresponding  Rayleigh-Schrocdinger  terms,  as  long  as  the  expansion  142 
converges  quickly. 

assumption  that  AE  is  small  relative  to  all  excitations  energies  in  the  "Q  space". 


(152) 
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In  the  normal  Mdller-Plesset  partitioning  138.  65.  69.  70,  89] 
H„"P  = fnP'l 

V"p  = i(p7l|«)pV«--0>HI?t)P,9 

E\,P  = 


In  general  E'UP  is  as  large  or  larger  than  the  total  energy.  However,  we  can 
introduce  a constant  shift 


c = = e'up 


redefining  II 0 and  1'  as 


H0  = H,"'’  + C 

V = vup  - c 


(153) 


(154) 


Using  this  choice,  reminiscent  of  the  Epstein-Nesbet  partitioning  [47,  88). 


£°  = Eur 

E'  = 0 


(155) 


Summarizing,  using  this  partitioning,  modified  Rayleigh-Schroedingcr 


(MRS) 


E?URS  = E?RS  + E'KS  = EHF 


(156) 


<157) 


El 


nn\abmmab\\ki) 

«»K‘»  +«!-«.-  <») 


+ (<,  + £>-(„-  £l)(<»  + £!-£«-  ‘it 


<,j||afc)(gfc||ij> 

(«,  + «»)"+' 


=J> 


(158) 


<159) 


(160) 


(161) 


(162) 


(163) 


1 (ullat)  (ai||crf)(c(f||0) 

s («. <.)"+'(<, + c,  - Cc  - ^r-*1 
1 <0||a6>(t/||.J>(o6|H/> 

8 (ti  + ‘s>m+1  (164) 

(ij\\abUak\\ic)(cb\\kj) 

(«.  + 1,  - n - n)"+1(«s  + n-n-  nr*' 

= y~l  fin  (E}lBW  + £.UBll')” 

The  lerm  of  lowest  exponent  in  the  Brillouin-Wigncr  expansions  (161-164) 
are  just  the  corresponding  Rayleigh-Schroedinger  expressions  (156-160). 

One  can  estimate  how  fast  the  polynomial  coefficients  decrease  by  noting  that 
d+ij-n-n  >2  6.  where 6 is  the  homo-lumo  gap.  leading  to  2gi  ~ fg1  - fjr- 
Hence  in  situations  where  the  homo-lumo  gap  is  not  too  small  (near  degeneracies), 
one  can  expect  a quick  convergence.  In  other  cases,  quasi-degenerate  theories 
should  be  used. 


Penurbation- Variational  Theory 

The  inclusion  of  just  a single  variational  parameter  in  front  of  the  first  order 
Rayleigh-Schroedinger  wavefunction  proved  to  give  very  good  results  even  in 
regions  of  the  potential  energy  surface  in  which  PT  does  not  converge. 


so 

The  advantages  of  perturbation-  variational  methods  are  their  far  greater 
speed  over  conventional  variational  methods  such  as  (matrix  diagonal ization) 
Configuration-Interaction  (Cl),  and  their  far  greater  stability  over  Perturbation 
Theories  for  larger  regions  of  the  Potential  Energy  Surface.  These  methods, 
however,  are  not  size  extensive  leading  to  conflicting  opinions  on  their  utility 
(19,  91-96). 

Using  the  approach  of  section  2,  one  can  write  the  first  order  Brillouin-Wigncr 
wavefunction  as 


|BIV,)  = ^(-£8IF)'I"  + 1)  065) 

\n)  = (E.-H')-°QVup\0) 

= (Rks)"QVUP\0)  (166) 

= ..S'  (°6||tj) 

l,; 

we  examine  the  following  approximation  to  the  ground  state  wavefuncuon 

|ti')  = C0|0)+  Y,  (167) 

where  the  coefficients,  {C,(.  are  variational  parameters,  and  NEXP  is  the 
number  of  the  “expanded  space  functions". 

This  expansion  is  a Krylov-like  expansion  (97).  In  the  limit  of  NEXP  = 
number  of  double-excitations  with  the  same  symmetry  as  the  ground  slate,  it  yields 
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CI-D.  However,  near-linear  dependencies  appear  well  before  the  true  theoretical 
linear  dependency  limit  is  approached. 

Near  dissociation,  where  usually  a few  specific  configuration  become  impor- 
tant. we  also  considered 


W = Cold)  + £ C.l»)  + X>«lw>  <168> 


where  P = |0)(0|  + g l£f)<jfl  and  |J?)  = cWffc |0)  is  limited  to  a subspace  of 
doubles  important  to  describe  the  bond  breaking.  This,  as  expected,  improves  the 
accuracy  and  the  numerical  stability  of  the  procedure,  but  at  the  cost  of  a modest 
amount  of  increased  computer  expense.  Below  we  prove  the  validity  of  using  a 
sequence  such  a s£  C,|u„)  = £ C„(4)"|i)  in  solving  the  Schroedinger  equation. 

We  wish  to  solve  the  eigenvalue  problem 


m=m 


B = B' 


(169) 


assuming  an  appropriate  |f»)  of  the  form  given  in  equations  167  or  168.  Let  us 
consider  the  operator  A.  and  its  non-degenerate  eigenvalues  and  eigenvectors 


i|«m)  = Qm|Om).  " 
A = A ' 


(170) 


We 


of  the  form 


K>  = (<*)». 


(171) 
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and  we  will  show  that  ihe  sel  (Inn)}?'  is  linear  independent;  i.e.  £ i?nK)  = 
0 =>!/„=  0,  Vn. 

To  prove  this,  we  resolve  |t>)  in  {|nm))f(  as  |t>)  - E^  Cml°m)  = 
jl  l«m){«m|v>. 

One  obtains 

0 = ^ hnM  = Yi 

= Y *>(*)"  Y. 

n=l  ^ ™=t  (172) 

= Y Z 1-<m{am)’\am) 

= ECml  l<im>X>..(Om)" 

E is  a number.  Thus  E CmKml«m)  = 0 =»  (.«.  = 

0,  Vm,  since  {|om)}J'f  is  a linear  independent. 

For  every  („^0ve  must  have  nm  = 0.  Let  us  assume  we  have  L < M 


X)ln(Oi)"  = 


(173) 


Thus  if  N < L =>  i)n  = 0,  Vn.  So  we  can  generate  a linear  independent  set 
for  the  N-  dimensional  space  spanned  by  { I t’.i ) } jV  = {(•d)"|0))i'  ■ 

In  this  space  we  want  to  solve  (B  - Bi  )|i>)  = 0.  Expanding. 

|4)  = X>.M  <l74> 

we  obtain  in  standard  fashion 

dtl(Bn  - 6 Sh)  = 0 (175) 


Bit  = (vi|B|t'i) 

fc-hW 


(176) 


In  our  present  case  \v)  = Q VM,,|0)  and  A = Rrs  = (£>  - Ho)~'Q-  50 
{|v„)}  = {|n)J  = {(B«s)"<?  V*,/'|0)}.  which  expands  the  double-excited  space 
of  same  symmetry  as  the  reference  function. 


AH  ihc  calculations  were  performed  using  INDO/1  Model  Hamiltonian  [81. 
82)  and  minimum  basis  set.  All  results  are  shown  in  atomic  units,  unless  explicitly 
said  otherwise. 

In  tables  3-7  , we  show  the  2nd  and  3rd-order  Brillouin-Wigner  results 
using  the  expansions  described  in  section  2.  for  water,  methane,  formaldehyde, 
acethylene,  and  benzene  near  the  experimental  equilibrium  geometries.  Nma»  is 
the  order  of  polynomial  in  the  right-hand  side  of  equation  162.  One  sees  that 
the  series  quickly  converges.  The  percent  error  is  with  respect  to  the  larger  Cl 
calculation  we  performed  . in  general  CI-Full. 

We  also,  can  see  the  importance  in  choosing  the  constant  C.  The  second 
column  shows  results  for  EgW  when  C = 0 , while  the  third  column  shows 
them  for  when  C = -j(W||W).  When  C = 0.  V = V'IF.  the  sequence  is 
not  convergent.  In  all  the  cases  examined,  the  zeroth  iteration,  the  RS  result,  is 
closer  to  the  CI-Full  result.  In  the  case  of  benzene,  we  do  not  have  the  Cl-Full 
result.  The  toughest  case  was  benzene,  which  took  a polynomial  of  7-lh  degree, 
however,  it  is  still  a easy  task. 
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Table  3:  H;0  — Brillouln-Wjgncr  Results 


h2o* 

Ehf= 

-18.120732 

AEci-d= 

-.038530 

(251) 

AEa-FulF 

-.039185 

Nm" 

£|,i 

(V  = P"p) 

E%w 

E\w 

E\w  + E%w 

0 

(RSPT) 

-.026250 

-.026250 

-.008457 

-.034707 

(11%) 

1 

-.214201 

-.025990 

-.008229 

-.034219 

2 

-1.42935 

-.025993 

-.008234 

-.034227 

(13%) 

a)  d|io=1.0130  A.  «hoh=104.5230°. 

Table  4:  CK,  — Brillouin-Wigner  Results 

ch4* 

Ehf= 

-9.260282 

AEci.d= 

-.042860 

(2%) 

AEci-Full= 

-.043689 

Nm« 

Elxv 

(V  = V-'IP) 

E2bw 

E%w 

0 

(RSPT) 

-.029022 

-.029022 

-.009682 

-.038701 

1 

-.171979 

-.028700 

-.009391 

-.038091 

2 

-.812863 

-.028703 

-.009397 

-.038100 

(13%) 

a)  dcH-1.0550  A,  «hch=109.47°. 
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Table  5:  H^CO  — Brillouin-Wigner  Results 


H2CO  ■ EHF=  AEa-D=  AEa-Fun= 

-24.441506  -.100381  -.106545 

(6%) 

Nmu 

Elw 

(V  = V"p) 

E\w 

ebw 

eBW  + ebw 

0 

(RSPT) 

-.080365 

-.080365 

-.015918 

-.096283 

(10%) 

1 

-1.30409 

-.077366 

-.014250 

-.091616 

2 

-10.9855 

-.077487 

-.014394 

-.091881 

3 

-21.5074 

-.077482 

-.014382 

-.091864 

(14%) 

a)  cicn=  1 -0990  A,  dco=l-203  A,  »hch=12I.75“. 


Table  6:  C2  H2  — Brillouin-Wigner  Results 


C2H2  * Ehf=  AEci-d=  AEci-foIi= 

-13.574992  -.103317  -.111199 

(7%) 

N"“ 

Eiw 

(V  = VMP) 

£flif 

£fltr 

E%w  + Esbk 

0 

(RSPT) 

-.089220 

-.089220 

-.012832 

-.102052 

(8%) 

1 

-.839939 

-.085128 

-.011383 

-.096511 

2 

-5.29892 

-.085324 

-.011517 

-.096841 

3 

-10.8770 

-.085314 

-.011505 

-.096819 

4 

-13.6370 

-.011506 

-.096820 

(13%) 

a)  dcc= 1.2024  A.  dcH=l-0750  A.  linear. 


Table  7:  Cs  H<s  — Brfflouln-V 


In  table  8 , we  can  see  that  the  variational  procedure  described  in  section 
3 soon  saturates.  With  the  exception  of  CH4.  all  other  cases  shows  instabilities 
when  the  lowest  and  largest  eigenvalues  of  the  overlap  matrix  S is  of  order 
ltr12.  Since  increasing  the  basis  set  must  monotonically  decrease  the  energy,  any 

of  5 expanded  functions  (NEXP=5)  seems  to  give  good  and  reliable  results  at  all 


from  the  theoretical  lie 


58 

computational  effort  is  the  same  as  3rd  order  Perturbation  Theory. 


Table  8:  Perturbation  - Vs 


1)  d|io=1.0130  A.  9hoh=1W  5230°. 
it)  doi-l-QSSO  A,  flHCH=109.47». 


d)  dcc=1.2024  A.  dcH=10750  A.  linear. 

e)  dcc-I-4000  A.  dcit= 1-0880  A. 
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Table  9:  HjO.  symmetric  sireich,  0hoh=1W-523O° 
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Table  11:  C2H2.  dot*  1-0750  A.  Unear 


Raylcigh-Schrocdinger  for  bond  lengih  and  force  constant.  Our  variational  method 

Figures  4-6  show  the  symmetrical  stretch  curve  for  water  at  a fixed  angle 
^HOH— 104.52°.  One  can  see,  from  figure  5,  that  around  equilibrium  all  approxi- 
mations are  fairly  good.  At  long  distances  (dno>  4.0  bohr)  Rayleigh-  Schrocdinger 
perturbation  theory  breaks  down  expected  (figure  4).  but  around  the  minimum  the 
Brillouin-  Wigner  results  closely  follow  die  Rayleigh-  Schroedinger  results.  Us- 


perturbation  theory  up  to  dHo=  3.5  bohr.  NEXP=5  (NEXP5)  follows  closely 
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results  obtained  for  the  most  computadonal  demanding  CID,  and  CI-Full.  Bet- 
ter agreement  is  obtained  when  singles  are  included  (NEXP5S).  since  the  bond 

NEXP5  is  stable  even  when  Perturbation  Theory  breaks  down,  again  as  expected. 

Figure  7 shows  the  breaking  of  the  double  bond  of  formaldehyde  using  the 
different  perturbation  theories,  second  order  Brillouin-  Wigner  and  third  Older 
Rayleigh-  Schroedinger  yields  the  same  force  constants  (table  10).  which  is  in 
better  agreement  with  the  full-CI  results.  In  figure  8.  we  sec  an  artificial  maximum 
which  appear  for  CID.  This  is  absent  in  NEXP5,  but  reappears  in  NEXP=5  + CID 

to  virtual  orbitals  (7-8)  (NEXP5.48).  It  should  be  noted  that  in  the  dissociation 
limit  (6.5  bohr).  CID  over  this  restricted  space  (CI-D.48)  and  NEXP5.48  present 
an  error  with  respect  to  CID  of  .004%,  while  NEXP5  presents  .291.  On  the  other 
hand,  near  the  equilibrium  geometry  (2.495  bohr)  CI-D_48  has  an  error  of  .291 
while  NEXP5  and  NEXP5.48  have  only  a .01%  error  with  respect  to  CID.  Hence 
NEXP5.48  seems  to  follow  closely  the  CID  curve  at  all  intemuclear  distances. 
The  inclusion  of  single  excitations,  figure  9.  removes  the  artificial  maximum  in 

Figure  10  and  1 1 show  the  breaking  of  acethylene's  triple  bond.  Perturbation 
Theory  breaks  down  very  early  (figure  7).  and  at  dcc>  3.0  bohr.  we  can  no  longer 
make  our  expansion  for  Brillouin-Wigner  perturbation  theory  converge.  Near  the 


Figure  10  and  lable 
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minimum  Brillouin-Wigncr  results  seem  beller  behaved  (see 
1 1).  Perhaps  this  is  due  to  the  fact  that  3rd-order  Brillouin- Wigner  results  are 
upper  bounds,  suggesting  an  easy  and  cheap  way  to  evaluate  the  reliability  of 
perturbation  theory  results.  However,  unlike  Rayleigh-Schroedinger  perturbation 
theory,  it  is  not  a size-extensive  theory. 

NEXP5  follows  CID  closely.  NEXP5.48  (NEXP5  + CID  over  a restricted 
space  spanned  by  double  excitation  from  occupied  orbital  (4-6)  to  virtual  orbitals 
(7-8) ) gives  very  good  results.  Again,  the  error  in  the  correlation  energy  expanded 
by  CID  is  0.1*  for  CI-D_38.  0.4*  for  NEXP5.  .01*  for  NEXP5_48  at  large 
intemuclear  distances  (4.5  a.u.);  and,  respectively.  0.3*.  0.05*.  and  0.02*  near 
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Figure  4:  H2O  — Perturbation  Theory 


Figure  5:  H2O  — Near  Equilibrium 


rgy(Hartrees)  Energy(Hartrees) 
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dHO(Bohr) 


dro(Bohr) 


Figure  7:  CH;0  — Penurbaiion  Theory 
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figure  9:  CHjO  — Single 
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Figure  II:  C;H: 


Summary  and  Conclusions 


Wc  have  described  a simple  expansion  of  lhe  Krylov  type  and  have  used  it  to 
derive  energy  expressions  for  Brillouin-Wigner  perturbation  iheory  from  Raylcigh- 
Schroendinger  perturbation  theory  as  well  as  define  a variation-perturbation  proce- 
dure using  a reduced  space  of  double  excitations,  thus  a "reduced  space  method". 

Brillouin-Wigner  perturbation  theory  calculations  take  very  Utile  additional 
time  over  the  Rayleigh-Schroendinger  perturbation  theory  calculations.  In  the 
examples  we  reported  here,  and  others  not  reported,  third  order  Rayleigh- 
Schroendinger  theory  recaptures  more  of  the  correlation  energy  than  docs 
Brillouin-Wigner  theory,  but  force  constants  and  bond  distances  are  better  re- 
produced using  second  order  Brillouin-Wigner  perturbation  theory. 

The  reduced  space  variational  procedure  NEXP5  follows  CID  closely  over  all 
regions  of  the  potential  energy  surface;  in  the  cases  examined  here,  the  difference 
in  the  correlation  energy  is  less  than  7%.  Variational  NEXP5_C1D  follows  CID 
closer  but  with  increased  computational  cost.  NEXP5S  iheory.  that  includes  single 
excitations,  seems  quite  accurate  when  compared  with  CI-Full,  and  it  is  this  quite 
efficient  theory  that  we  now  examine  for  electronic  spectroscopy. 


CHAPTER  4 

FREQUENCY  INDEPENDENT  RESPONSE  PROPERTIES 


Introduction 

We  investigate  nature  through  the  response  of  a system  to  some  external 
factor.  One  of  the  most  common  probes  is  a static  electric  field  (E).  Let  us 
investigate  165)  the  response  of  a state  |d-)  corresponding  to  H to  an  external 
time-independent  one-electron  perturbation  AHj. 

H - H + AHi  (177) 

E(  A)  = $\H  + A//j|ti)  = Eo  + \Ei  + ^A  2E3  + ^A3£3  + . . . (178) 

The  terms  that  are  nonlinear  in  A arise  because  the  state  wavefunction  Id1)  depends 
on  A.  or  has  "responded"  to  the  perturbation.  When,  for  example,  AHi  represents 
a static  electric  field  (AHi  = E.r).  Ei  yields  the  permanent  electric  dipole  moment 
(p)  of  the  unperturbed  state  |i>),  £2  gives  the  state's  polarizability  (o),  and  E3.E<. 
etc.  yield  higher  hyperpolarizabilitics  (0,  7, etc.). 


The  presence  of  an  electric  field  can  dramatically  change  the  molecular  orbitals 
of  the  system.  The  inclusion  of  the  single  excitation  manifold  does  not  change  the 
zero-order  energy,  but  its  presence  allows  us  to  account  for  the  orbital  relaxation 
due  to  the  field.  However,  the  calculation  of  properties  is  often  more  sensitive 
to  the  quality  of  the  wave  function  than  the  energy,  so  the  use  of  a correlated 
zero-order  energy  seems  appropriated. 

We  propose  two  ways  to  accomplish  the  inclusion  of  relaxation  and  correlation 
in  the  calculation  of  dipole  moment  and  polarizability.  The  first  is  reminiscent  of 
a CC-S  applied  to  a correlated  reference  function.  The  second  can  be  viewed  as 
a linearized  CC-S.  or  more  appropriately  a Tamm-Dancoff  Approximation  (TOA) 
or  Cl-S  on  a correlated  reference  function. 

In  section  2.  we  present  a short  review  on  CC  response  properties.  In  section 
3 . we  introduce  our  modified  CC-S  and  CI-S  response  properties.  Finally,  we 
summarize  out  results  and  present  our  conclusions  in  section  4. 

The  Coupled-Cluster  Response  Properties  1 65.  98] 

Assuming  that  the  Hamiltonian  can  be  written  as 

tf(  A)  = « + XH,  (179) 


and  the  wave  operator 


n(A)  = exp|r(A)) 


(180) 


(pressed 


|0(A))  = fi(A)IO)  (181) 

exp[-r(A)l//(  A)  exp[7-(A))|0)  = £(A)|0>  (182) 

Neglecting  the  variations  in  the  orbitals  when  the  field  is  applied,  the  equations 
for  the  energy  and  cluster  amplitudes  become 

(0|exP|-7-<A)l«(A)exP[r<A))|0)  = £(A)  (183) 

(•jn«P(-rW|W(A)exp[r(A)]|0)  = 0 (184) 

The  above  Hamiltonian  can  be  expanded  in  powers  of  the  field  A once  the 
duster  operator  T( A)  and  H( A)  are  so  expanded: 

T(A)  = T°  + ATl1»  + AJrl2l  + ...  (185) 

£(A)  = £°  + A£(i,  + A!£,j»+...  (186) 


The  resultant  first-  and  second-order  equations 


£<‘l  = (0|e*p[-T°){W,  + [£,T<1>])exp|r°]|0)  (187) 

0 = <Sjpexp|-T0H«i  + |£,r1,|])«p[7J>]|0)  (188) 

£(J)  = (0|exp(-r0]{|//,,T<1,]+5ll«.r|1,],r"l]+[W,T<2»]}exp[7-°]|0)  (189) 
o = (5S:::|exp[-r°]{[£,.r<l>]+|[(£,r',)],r<I>i+|£,r(I>]}exp[T0)|o)  (iso) 
Dipole  Moment  and  Polarizabilit)1 

Within  the  Cl  framework,  the  dipole  moment  is  just  an  expectation  value  of 
the  ground  state,  hence  the  dipole  moment  calculation  is  a direct  measure  of  how 
"good"  the  wavefunction  is.  On  the  other  hand,  one  can  view  the  expectation 
value  of  the  one  electron  operator  as  a sum  of  the  overlap  between  the  Cl 
wavefunction  and  a function  which  includes  single  excitations  and  de-excitations 
with  respect  to  the  Cl  ground  state.  Coupled  Hartree-Fock  (CHF)  or  Random 
Phase  Approximation  (RPA)  needs  only  single  and  double  excitations  to  describe 
dipole  moment  and  polarizability. 
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Oddershede  and  Svendsen  [99]  have  shown  that  the  polarizability  tensor  is 
well  described  within  the  CHF  approximation. Ohm  et  al  [100J  have  observed 
that  the  extension  of  RPA  and  TDA  to  a correlated  reference  function  results  in 
similar  results  for  transition  energies  and  transition  moments.  Thus,  we  propose 
to  use  a wavefunction  that  is  a combination  of  single  excitations  on  a correlated 
reference  function. 

The  first  method  is  similar  to  a Coupled  Cluster-Singles  (CC-S)  applied  to 
a first  order  perturbation  theory  reference  function,  or  a perturbation-variational 
function. 

The  second  method  can  be  viewed  as  a linearized  CC-S.  where  only  the  terms 
depending  linearly  of  the  cluster  amplitudes  are  retained.  This  method  can  also 
be  interpreted  as  a generalized  Tamm-Dancoff  Approximation,  or  Cl-S.  where  the 
reference  function  is  corrected  to  first  order  in  perturbation  theory,  or  is  obtained 
through  a perturbation-variational  procedure. 

Both  methods  are  solved  requiring  that  the  Schroedinger  equation  be  obeyed 
in  the  subspace  spanned  by  the  Hartree-Fock  wavefunctions  and  its  single  excited 
configurations. 


Explicitly. 


Modified  CC-S 


Assuming  Hi 

IW  = |0)  + |1) 

|1)  = flV|0> 

Efo  = (Oil'll)  = (0|V£V|0) 

Q re  = £.  RSPT 
(£-//,)  l £ = £ BWPT 
Q=  1 — |0)<0|  — |f>(*|  091) 

= fnP'l 

V = (pgllrs)pVsr  - (pk\\qk)p'q 

W')  = «r,W;n  = i?a'i 
l!;D  = «'i|i> 

Response  Properties 

= £.(p|r)g)  p\  is  an  one-elecuon  operator. 

fl(A)  = H + A Hi 

*p|-ri(A)]ff(A)exp[r,(A))|W  = £(A)|t!'.) 

(192) 

7|(A)  = 7f1  + AT}1)  +AJr,(J)  + ... 

£(A)  = £*“'  + A £>'»  + A2  £'=)  + ... 

1. 

£<“>  = (0|e-7t',We75°V.)  (193) 

(Si^W^o 


(194) 


EV»  = Ear  + £^  + i“(0)/,a  + ilf  if  (195) 


Epr  ‘f  = A/  + (Jl»|l>+ 
if  JU  - if  - if  if  /,.  + if  <M||ja)+ 

<f{(5|//|fi  1)  + ‘f  (ifclloj)  - if  <«IM  - if  if  (»*||«e)} 

(196) 


£*‘1  = (0|e-I,O,|/fi  + (197) 


($|<r7W{/f1  + [//,r«]  }eT<0>| *,)  = 0 (198) 


£<‘>  = (0|£,|0)  + (0|if,  |“)  if + 

/i.‘f + WII«&)‘f‘f 


(199) 
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£W  (*'*>  = -£(»(* <»> 


(;|W,|0)  + (J|Wl|l)+ 

(5i//,i?)  <;,0,+ 


(200) 


£I3)  = (0|e-I,'’|[Wi,T"))  + i[[H,.7-'"],7-(1,l  + ttf.TWlJe^V) 
(201) 


(5lt-I<,,|[Hl,r»>]  + i|[//,,rl"],r«1>)  + [H,r,!,]}tI<0V.)  = o (202) 


£lIl  = (0|H,|?)(Jll|  + /,.  tfh 

(203) 
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£<°)  ,‘<2>  = _£(2)  j*<°>  _ fid)  («‘>+ 
(‘1/7,1?)  t“<1)+ 

[({|H|?)  + (}IwI?;i»('<2,+ 

(,k|W|?»)  " + (jl^l?^)  <“t0W>+ 

<;1!,«5i«i?f>  «f’ + (*i«i?«  ) 'fV0’) 


(204) 


The  appearance  of  lerms  connecting  singles  with  doubles  and  triples  makes 
this  approximated  method  not  very  attractive.  The  linear  terms  in  the  cluster 
amplitudes,  on  the  other  hand,  have  the  signature  that  we  are  looking  for.  So. 
we  introduce  the  Modified  CI-S. 


Modified  Cl-S 


Similarly  to  the  Modified  CC-S.  we  define 
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W = |0>  + 11) 

|1)  = «V|0) 

Err  = <0|V|1>  = <0|V^V]0) 

» Q fe  = Eo  RSPT 
(£-fl0)  U = £<°>  BWPT 
Q=i-|0)(0|-|?){?| 

>h  = Inv'l 

V = (p9||rs)pVsr  - (P*ll**)Pf« 
W = (1  + C?|)W 
Cl  = C'fl'i 
|J;1>  = a*«|l> 


(205) 


H(  A)  = H + XH, 

m\)(l  + C,(A))|*„)  = £(A)(1  + C,(A))W'o) 
C|(A)  = C|0»  + AC,1)  + AJC(J)  + 

£(A)  = £*“'  + A £">  + A!  Em  + . . . 


(206) 


78 

Collecting  all  terms  of  same  order,  one  obtains: 

A":  [/f-£(°>)[l  + C1(o)]|V.)  = 0 

A1  : [H-  E^lcj'Vo)  = -[£,  - £">]|1  + C,<0)]|tM 

A2  : | H-  E|0ljC,,2V.)  = -I«t  - £(1))  C}1  V.)  + £|2,l'  + Cf’lW 

A3  : \H  — £,0|]C[3)|Vo)  = -[«l  - E1'1]  CjJV.)+ 

E(2|c!‘V-.)  + £(3||1  + cf’lW.) 

(207) 

One  can  obtain  the  equations  for  (C|"'(  and  {£*">)  by  projecting  in  the 
space  spanned  by  |0)  and  {|°)}.  Explicitly. 

£<">  = Ehf  + Epr  + (208) 

£(0(o»(0|  ={»|W|o)(=  /,  ) + (‘|E|1)+ 

’ (209) 

C‘|0){(‘|E|“)  + (*|W|;;1)} 

El"  = (0|W,|0)  + C'(0|(0|E,|f)  + C;"1/..  (210) 

£(a>c‘0)  = _ B(DC*I0)+ 

<;i//,io>+(}i«1|i>+ 

(211) 


Em  _ c,,|l|{o|Hi|“)  + c’mf,Q 


(212) 


£(0)C*P)  = _ £(‘)c‘(l1  - £<2>c‘(0|+ 

C"(1|(5ltf,lf)+  (213) 

c,°(J,{(*|W|“)  + (>|W|?;i)( 

E(3)  m C1,(2'(0|W1|“)  + (214) 

E,0lc*m  = - eCIcJ*21  - El'lc*1"  - £,3)cJ(0)+ 

C'<2)(‘|H||J)+  (215) 

C*(3,«;i«l?)  + (*|ff|?;l)} 

Specifically,  one  obtain  for  the  dipole  moment  in  the  direction  t'l 

= <0|x,|0>  + (Olx.it)  cfm  (216) 

and  for  the  component  i j of  the  polarizability 

0,2  = + (2,7> 
and  for  the  component  i j.k  of  the  hyperpolarizability 

02^4(2)  a2£,d(2)  n2£*4(2) 

a,‘l  = dE,‘dE,^Xt^  + d£18£t*0|ljl')  + a£J8£t*0|r'1^  <218) 

where  E, . . £*  are  the  fields  in  the  directions  t'i  ,cj  and  ea,  respectively. 


(219) 


£<°>c:«01  = mm + mmi) + ami-,  «j  c*w  (220) 

E“'^Kr 

tort  <5""+ 

#hO^+aws^+  » 


in  A3. 


The  CI-SD,  CI-SDT,  CI-F  and  the  reduced  expanded  space  results  (NEXP) 


were  obtained  by  interpolating  seven  points  to  a third  order  polynomial 

W{E)  = Wo  + fE  + ^E,  + ^Et  (223) 

The  electric  field  strengths  were  0,  ±.001  a.u.,  ±.002  a.u„  ±.003  a.u.  for  the 
principal  axis  of  the  molecule;  and  0.  ±.01  a.u.,  ±.02  a.u„  ±.03  a.u.  for  the  other 
directions.  The  mean  square  deviation  was  ± 10"’. 

In  tables  12-15  we  present  the  results  for  the  Modified  CI-S  (MCI-S)  obtained 
from  equations  219-221.  Modified  CI-S  without  triples  (MCI-S_WT),  i.e.  without 
the  terms  in  equations  219-221.  We  also  present  the  NEXP  results 

obtained  from  applying  the  variational  principle  to  a wave  function  of  the  form 
W = |0)  + C1|1)  + C»|J)  (224) 

where  |1)  = /?«sl'|0>  and  |“)  = ^(oji0  + o^t«)|0).  The  polarizabilities  arc 
obtaining  by  interpolation  as  described  above  for  Cl. 

We  compare  these  results  with  Harltee-Fock.  RPA.  Cl-S,  CI-SD.  CI-SDT  and 
CI-F  for  water,  carbon  monoxide,  formaldehyde,  and  acethylene. 

By  Hartree-Foclc  we  refer  to  a sum  over  states  formula  with  Moller-Plessct 
energy  denominators.  Ci-S  is  also  a sum  over  states  approach,  with  the  sum 
over  diagonal  states.  Both  might  be  considered  as  approximations  to  the  RPA 
theory,  which  is  the  full  response  of  the  Hartree-Fock  wavefunction  to  an  ex- 
ternal perturbation.  CI-SD,  CI-SDT.  CI-F  and  NEXP  are  the  full  response  of 


82 

Ihesc  wavefunciions  lo  the  external  field.  The  MCI  results  discussed  below  arc 
appoximalions  to  the  Cl-F  results,  as  are.  in  a real  sense,  all  the  results  repotted. 

From  the  polarizabilities  results  we  see  that  Cl-SD  yields  better  results 
than  CI-SDT,  except  in  field  directions  which  perpendicular  to  the  plane  of  the 
molecule.  This  is  also  followed  by  our  MCI  results.  Although  the  overall  error 
between  MCI-S  and  MCI-S.WT  are  about  the  same  ±10%.  the  MCI-S_WT  seems 
to  describe  all  directions  better  than  MCI-S.  The  NEXP  results  closely  follows 
Cl-SD  in  all  cases. 

RPA  also  performs  remarkably  well.  This  is  due  the  very  balanced  way  in 
which  single  and  double  excitations  are  included.  However,  for  acethylene.  the 
errors  for  the  different  directions  of  the  fields  (table  15)  is  15%  for  the  parallel 
direction  and  only  2%  in  the  perpendicular  directions,  while  MCI*S_WT  yields  a 


From  the  dipole  moment  results.  Cl-SD  also  seems  to  outperform  Cl-SDT 
but  in  a less  dramatic  fashion.  This  is  followed  by  our  MCI-S  and  MC1-S_WT. 
The  average  error  for  Cl-SD  is  4%.  for  CI-SDT  is  5%.  for  MCI-S  is  10%  and 
MCI-S.WT  is  15%. 

In  conclusion,  it  seems  to  be  more  accurate  to  use  RPA  or  Cl-SD  than  CI-SDT 


1-S_WT  is  a very 
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yields  polarizabilities  within  ±10%  of  the  CI-F  results  and  it  is  an  N4  method. 

are  done.  In  this  framework  the  NEXP  method  is  competitive  with  CI-SD,  since 
this  method  goes  with  N5,  while  CI-SD  dependence  is  N6. 
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Table  13:  Polarizability  for  Carbon  monoxide. 
dco=l-12819  A . Ehf=-  22.619953  Harirees. 


Method 

li  (Debyes) 

Q„  (A3) 

a,  (A3) 

-.0101 

-1.9071 

-.4358 

(98%) 

(33%) 

(47%) 

-2.0491 

-.9467 

(43%) 

(15%) 

-1.4455 

-.8318 

(1%) 

(0.8%) 

-.5246 

-1.4713 

-.8263 

(11%) 

(3%) 

(.1%) 

-.5142 

-1.3866 

-.8269 

(13%) 

(3%) 

(.2%) 

-.5130 

-1.4289 

-.8197 

(13%) 

(0%) 

(.6%) 

-.4679 

-1.7436 

-.8263 

(20% ) 

(22%) 

(0.1%) 

-.3887 

-1.4629 

-.7255 

(34%) 

(2%) 

(12%) 

Cl-F 

-.5892 

• 1.4289 

-.8251 

t 
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CHAPTER  5 

REDUCED-EXPANDED  SPACE  METHOD 
APPLED  TO  ELECTRONIC  SPECTROSCOPY 


Iniroduclion 

Forcsman,  Head-Gordon,  Pople  and  Frisch  1102]  divide  eleclronic  structure 
theories  for  excited-states  chemistry  into  two  types: 

1.  Those  methods  which  concentrate  on  the  physics  of  the  traasition.  e.g.  Ran- 
dom Phase  Approximation  (57.  65,  103.  104],  and  other  Green  Function 
methods  (57,  65,  105,  106], 

2.  Those  methods  which  concentrate  on  the  physics  of  the  state,  e.g.  SCF,  Cl, 
MCSCF,  GVB,  SAC-CL 

They  also  suggested  a general  criteria  for  an  acceptable  theoretical  model  for 


defined,  applicable  to  systems  of  interest,  size-consistent,  and  variational.  Turning 

framework  should  lead  to  accurate  representations  of  several  electronic  states,  and 
these  should  have  wave  functions  which  are  directly  comparable.  This  qualification 
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is  necessary  in  order  I o calculate  transition  properties  among  the  various  possible 
states.  It  also  ensures  the  ability  to  quantify  the  differences  between  the  structure 
of  the  excited  and  ground  slates... 

In  this  chapter  we  investigate  an  approximate  scheme  that  tries  to  fulfill  all  of 
Pople's  criteria  with  the  exception  of  size-consistency,  while  trying  to  incorporate 
the  physics  of  the  state  and  the  transition.  Also,  we  keep  strong  ties  with  TDA, 
which  allows  for  an  easy  interpretation  of  the  spectra. 

In  section  2.  we  describe  our  method  for  excited  states.  In  section  3.  we 
present  results  for  some  selected  systems:  water,  formaldehyde,  and  acetylene. 
Finally,  in  section  4,  we  present  a summary  and  our  conclusions. 

Reduced-Expandcd  Space  Method  applied  to  excited  states 

The  lowest  level  of  theory  to  yield  excitation  energies,  transition  moments, 
and  wave  functions  is  single  excitation  configuration  interaction  (CI-S).  or  the 
Tamm-Dancoff  approximation  (TDA).  This  approximation  consists  in  using  the 
orbitals  of  a HF  state  in  an  otdinary  Cl  procedure  and  solving  for  higher  roots. 
Excited  determinants  are  produced  by  replacing  one  occupied  orbital  by  a virtual 
orbital.  Considering  the  simplicity  of  this  approximation,  the  results  are  quite 
remarkable.  One  generally  obtains  the  right  order  of  the  low  lying  transitions, 
and  fairly  good  transition  energies,  however,  the  transition  moments  are  often 
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exaggerated.  One  can  greatly  improve  the  transition  moments  by  the  inclusion 
of  double  excitation  to  the  manifold  (CI-SD);  however  this  usually  makes  the 

with  the  double  excited  slate  manifold  which  is  not  balanced  with  respect  to  the 
single  excited  states.  Sometimes  the  inclusion  of  double  excitations  even  reverse 
the  correct  order  of  the  slates.  The  inclusion  of  triple  excitations,  specially  those 
which  double  excitations  with  respect  to  the  single  excited  stale,  corrects  this  by 
lowering  the  energies  of  the  single  excited  states.  Schematically. 


RPA  |103]  usually  gives  good  transition  energies  and  transition  moments. 
This  accomplished  by  a careful  inclusion  of  some  double  excitations.  However, 
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one  looses  the  description  of  stales.  In  fact,  one  can  inleiprel  the  RPA  as  if  each 
excited  state  had  a ground  state  especially  derived  to  insure  the  calculation  of 
good  oscillator  strengths. 

So.  one  concludes  that  in  the  chemistry  of  excited  states,  it  is  desirable  to 
include  correlation  in  older  to  obtain  reliable  transition  moments.  However,  it 
must  be  included  in  a balanced  fashion;  i.e.  all  states  should  be  correlated  at 
the  "same  level"  to  avoid  distortions,  as  occur,  for  example,  in  such  procedures 
as  CI-SD. 

In  many  cases,  the  correlation  present  in  the  ground  state  and  in  the  excited 
slates  are  very  similar.  In  such  cases  one  can  apply  single  excitations  to  a reference 
function  optimized  for  the  ground  state  and  obtain  a good  approximation  of  a 
correlated  excited  state.  One  can  say  that  the  correlation  is  portable.  This  is 
the  idea  behind  SAC-C1  [1-18]  which  has  proven  quite  successful.  And  it  also 
confirms  the  observation  of  Ohm  and  co-workers  [ 100]  that  RPA  and  TDA  on  an 
antisymmetric  geminal  power  reference  function  yields  similar  results. 

We  propose,  for  a closed-shell  Haruee-Fock  reference  function  (|0>),  a mul- 
tielectron basis  set  of  the  form 

{|0).{l?)}i{l!;)}i  {|zV)J,  {|°;  A')})  (225) 

explicitly 

|“)  = a'i|0)  are  the  single  excited  configurations. 
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||f)  = c'i/'R- 10)  are  some  specific  double  exciled  configurations, 

|AT)  = _'a  jw <?V |0) ; N = l.A’mu  are  the  functions  which  represents  the 

double  exciled  configurations  manifold  not  explicitly  included. 

The  triple  excitations  are  included  in  a similar  fashion  to  the  double  excita- 
tions. We  express  this  in  two  different  ways:  the  first. 

|f;  N)  = (£  ‘^jtfQa'igVIO ):  N = l.AW  026) 

has  a denominator  which  is  the  sum  of  the  orbital  energies  of  three  holes  minus 
the  sum  of  the  orbital  energies  of  three  particles;  while  the  second. 

If;  N)  = o'iQ ! fQV|0);  N = 1 ,Nma,  (227) 

is  a single  excitation  on  the  "average"  double  exciled  functions 

0 = i - (|0)(0|  - |*){fl  - IhXhI’)  (22gj 

P + Q = i 

H = H„  + V (229) 


= /«pV.  /m  = (pl  - ^'s72k)  + (p*!!?*) 


(230) 
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V = j(w|MpV*r  - Ml?*)?’?  <23l> 

(P?llrJ)  = J ¥’I(*l)V>l(*»)^i|~V^(xlWx*)dxidxi  <232) 


|0);  (0|tf|0)  = E,if  (233) 


E.  = (0|H.|0)  = ftk  (234) 


(235) 


Hence  one  has  to  solve  a secular  problem  of  the  form 


(236) 


(237) 


H'  = H-{0\H\0) 


E,„,  = E - Eh  f 


(238) 


W-<*hW 


(239) 


Hij  = (MH'Wi)  (240) 


Explicit  expressions  for  ihe  elements  of  matrix 

If  one  neglects  the  terms  in  parenthesis  in  the  average  triples  denominator, 
one  gets  the  second  formulation  for  the  average  triples. 

Overlap,  A,,  = (i',|0,) 


(242) 


(243) 


(0||f>  = ' 


(244) 
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(0|  A')  = 0 (245) 

(0|“;A')  = 0 (246) 

(*l?)  = (247) 

(}ltf>  = 0 <248) 

(*|.V)  = 0 (249) 

(*|f;A')  = 0 (250) 

(flLlg)  = j(«/«V  - (251) 

(*  |JV>  = 0 (252) 


(&„I?;W  = o 


(253) 


[M\N)  = 


4 (£,  4 Cj  - £.  - £»)"+■' 


(254) 


V)  = 

Kit.,  mcMcdm 

4 Ui  + £|  - e,  - u + (£|  - £«))"+•'' 

((.-.||c4)M|(-j) 

2 (£,  + £4  - £C  - U + (4j  - £.))"(£*  + 4>-«c  - «£  + (£.  - £»))•' 

6^ (t(||CT)(ct||K) 

2 (e4  + £,-£,-£„  + (£,  - f4))W(£4  +£!-£<-££+  (fi  - £.))-'’ 
(ti||cc)(c6||tj) 

(£,  + £t  - £C  - e.  + (£j  - £*))*'(«£  + £4  - e,  - £4  + (£,  - £,  ))A 


Hamiltonian.  //’  = H — (Q|//|Q) 


<0|//'|0>  = 0 


(257) 


<0|«T)  = /..  = 0 


(258) 


<0|ff'IS?>  = (H||cd> 


(259) 


(26!) 


(}|#T)  = t.j/i.  ~ + (>t||oj>  (262) 


(}l«'ISf)  = 

+ + «,i(W|M)  (263) 

-6l'{kl\\}<l)-Sumcj) 


(‘l«W  = 

1 (HH^)MHv) 

2 (£4  4 e,  - ec  - Cj)x 

1 (kl\\cj)(d>\\kt) 

” 2 (£*  + £|  - Ec  - £»)* 


(264) 
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d’Jtf'lif)  = 

(6/effi  ~ 6/if ic  + 6,i!f,  - 6tcffi)(6 mtKl  ~ 

- (6ml f I*  - 6mifi„  + 6,iftm  - 6„tfim)(6f,6,j  - 6fj6,')+ 

(6mi6„i  - + (6/,6,j  - 8/j<tc)(wllmn)+ 

i/i(I-*mlr((jl|CT>  + + 4i,j(ij||«n)]+ 

Vl+4mt('/l|cn)  - + 6M\\cm)  - 6ni(lf\\cm)}+ 

4/c(+4mi('sll*0  - “ 4|.*('»II<M]+ 

6tcl-6mtmd«)  + - M*/I|rfm)  + 4„t(//||rfm>] 


(LWW  = 

1 (f g\\cd)(cd\\mi 


1 (*t||mn)(/g||*-f) 

2 (£t  + £(-£/-  £,)  V 
(*j||£m)(c/|H-n) 

(£t  + £„  - £C  - £/)  V 
(kf(\cm}(c9\\kn) 

' (si  + £„-£,-  C,)K 


(267) 


(e*  + £n 


<&,m*sA'>= 

(■g||ac)(/c||mn) </i||ea)(eg||mn) 


-£,  + (£,-£.))•'  (£t  + £„ 


1 2 (£m+£4-£e-£j  + (£j-£.))J' 

M./  (*j|M{efl|*m> 

2 (£*  + £„-£«-£</  + (£,-£„))•' 
(*/||crf}(cd||in) 

2 (£m  + £»-£c-££  + (£i  -£«))' 

' 2 (Et+£„-fc -££  + (£, -£«.■'•' 

(268) 

m«n){fg\\mk) 

(KllmaX/sIM 

(£»  + ££-£/-£,  + (£.-£.))•'■ 

(£t  +£„-£/-  £9  + (£i  - £«))A 

, 1 (W||»n«>(/j||W> 

, 1 (t/ll«»X/j||H> 

‘"2  (ft  + £,-£/-£„  + (£,’  - Cc))K 

■”2  (£t  + £,-£/-£,+  (£,  - £«))■' 

(i*||mc)(/c||ni-}  , 

+ £*  - £/  - £c  + (£.  - £«))* 

(e» +£»-£/-£«  + (£( -£»))•'' 

(it||ne)(jc||m«)  , 

(fm  + £4  — £j  — £r  + (£i  - £b))A 

K.t, » (W||cm)(c/||W) 

(H||cnKc/||*()  &,/*■  {kl\\cm)(cg\\kl) 

(Ei  + £|- £/ - Ec  + (£i  - E«))A  2 (ct  + £|  - £e  - £,  + (£i  - £«))' 

(JV'|ff'|JV)  = -{A’’  - 1|AT)(=  <JV'|Ar  — 1))+ 

i mmusMicdm  , 

8 (£i  + £(  - £/  - f9r  (£i  + £J  - £c  - Erf)* 

1 (mn|M>(M||mn)(cd||t()  | (274) 

8 (Em  + E»  - Ec  - £j)A  (£t  + £(-£<-  ej)A 

(£m  + £/-£/-  Erf)  (ft  + £|  - Ec  - £j)A 


(A"|W'|f;  A’)  = 

+ (■t||/c>(/t||ad)(erf|H-f) . 

(£,  + £(-£/-  £c)ff’(£t  + El  - Ec  - Erf  + <£■  - £.))'' 

+ i (■t||g/)(/(|MM|*Q . 

7 (£,  +£!-£/-  Ea)K  (£t  + £|  - £<  - Erf  + (£i  - £«))’' 

+ 1 (t/||/4)(/.|l«.Kcrf||t)) . 

2 (££  + £(-£/-  Erf)''  (£»  + S|  - Ec  - Erf  + (£.  - «•))* 

_l mfli)(f,\\cj)(cdm r 

4 (Et  + E|  — Eo  — £/)  (Et  + El  “ Ef  - Ed  + (£t  “ Erf))' 

(mk\\<,c)(,l\\md)(cJ\\kl) 

(Em  + Et  - Ec  - Ea)A"(c*  + E,  - Et  - £J  + (f,  - £„))' 

_ 1 (MM)MI°™)("W r 

2 (Em  + £*  - Ec  - ej)K\ck  + E(  - £c  - Cd  + (£.  - £.))' 

_1 (mi||ai)(t<||md)(crf||t<> ? 

2 (fit  + Em  - Ec  - *•)*’(«*  + E|  - Ec  - Erf  + (£.  - Erf))' 

+ i (mi||cdH*<||n»«)(crf||*/) 

4 (Em  + S.  - Ec  - Erf)  ’(+£»  + £|  - Ec  - Erf  + (e.  - Erf)) 
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(;i  N'\H'\°-,N)  = [-(*;  A'|?;  {N  — 1))+] 

{=  -({iW'lf'.tA'-  1)>  -Ui  -£.)(};  JV'|;;AT)}+ 

where  ihe  lerms  in  braces  are  for  the  |* ; A')  = o,i/?v  V|0)  Iriples  and  the 
terms  in  brackets  are  for  the  |";A’)  = flAViV|0),  above. 


4ij4i  (i mf9){f9\M(cmi) r+ 

“ ® (f*  + £(-£/-£,+  (£;  - £i))A  (£1  + £|  - £c  - ££  + (£i  - £«))* 

1 (mn||cd)(kf||mn)(crf|H-i) 

* (£.«  + £«  - £c  - £a  + (£j  - £»))i#'(£4  + £(-£,-££+  U.  - £«))"' 

(m)||/rf)M|mc>M|*Q t)+ 

(£,„  + £(  - £/  - £J  + (£j  - Ej))A‘(£i  4 £|  - £f  - fj  + (f,  - £«))  ' 


4 (! (Wg)(/gMMI*>) r+ 

4 (=-.  + £4  -£/-£,  + {Cj  ~ £»))A'(£j  + £4  - £c  - £J  + <£|  ~ £.))"' 

1  (md)('f\\cj)(cdm _+ 

2 (£4  + £/  - U - £/  + Uj  - £4))A  (£4  + £)  - £c  - Ej  + (Si  - £«))* 

MIAQ{*/lfrQMI*>) r+ 

(£,  + £|  - £/  - £j  + (£,  - £»))A’{£>  + £4  - £r  - £J  + (fi  - £«))"' 

mmmcMcdm r+ 

(£(  + E|  - £/  - £J  + (fj  - £fr))  (£4  + £t  - £c  - ££  + (£i  - £. ))* 

mcd)(dmHfd\\lj) 

(£4  + E|  - C,  - £a  + (£,  - £4»A,'(£j  + £|  - £/  - £J  + <£.'  - £«))'' 

1 <*m  | |crf)  {<■  | |m j ) («/|  |*f ) 

2 (£4  + £,„  -£,-££+  (f,  - £s))V(£4  + £|  - £ C - £J  + U.  - £«))* 


1 (tf||cd)(im||t()(crf||jm) 

4 (ej  + £i - £c - ej  + (e,  - £i))iv'(ej  + cm-cc-ej  + (e<  - e»))a 

1 (im||erf)(H||jm)(crf||<r/) _} 

4 (£i  + e„  - ec  - Ci  4 (£,  - £j))A'(£i  + £j  - £t  - U + (£i  - e«))A 


4 (fm  4 E„  - £0  - f.  4 (e>  - £i))A’(£*  4 £1  - £ 
1 (Ur||m<.)M|M(oq|«) 

— Et4(s*— £<))* 

2 (£|4£m- 

c - Ej  4 (£,  -tl))A'(£«4£|-£ 

- ££  4 (e.  - £«))a 

(kd\\cm)(ml\\d<i){cb\\kl 

«-££  + (£/-  £»))*"(£*  + £1 " 

-£C  4 (£.-£.)/ 

(m(||0rf)Ui.||£m)WI|W 

(£m  4 £|  - 

a -£J4(£,  -£t))A'(££4£|- 

-££4(£j-£«))A 

(kl\\cd)(cm\\ka)(bd\\ml) 

~ £J  4 (£,  - £j))A'(£m  4 £|  - 

- U 4 (e,  - £.))A 

1 

<l-/||c/>(fc/||£.<#><C</||A-/> 

2 (££  + £,-£ 

- £/  4 (£j  - £i))A  (st  4 £/  - 

- £J  4 (£,  -£. ))" 

-£/4(E,-£l))A'(Et4£|-£ 

-£j  4 (£(-£.))" 

1 

4 (£*  + £(  - Cc  - u 4 (tj  - El)|A"  (ei  4 £|  - El  - £/  4 <£i  - «■„ )) r 


<»*ll/c)(6/M<to)  <erf»fc j> _+ 

(e,  4 £4  - £c  - £/  + (£j  - £i))  (£j  4 £i  - £c  - U 4 (£.  -£«))* 

+ 1 <ki\\cd)(cJ\\Ia)(bf\\jk) 

2 (£i  4 E*  -£«-££  4 (£,  - £»))A‘(e>  4 £1  -£»-£/  4 (£,  - £.))* 

1 maf)(mcJ)(cd\\kj) 

2 (£,  4 £t  -£.-£/  4 (£,  - £»))"'(£>  4 Ct  - £c 


- EJ  4 (£i  -£.))' 
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(imlNX/.lbmMcHItl) 

(f*  + £m  - £«  - «e  + (fj  - £»))"'(£*  + £/  - £J  - 
1 {im\\ac)(kl\MW\kl) 

-£.  + («-*.))" 

2 (e,  + £ro  - £,  - £C  + (e,  - £»))A"(£J  + £|  - £| 
1 (H||c«){.m||t/>(ei||;m) 

l-£c + (£,-£.))•'' 

2 (ej  + £,  - £„  - ec  + (t,  - £»))*'(£,-  + £m  - £. 

! 1 {H|M(«i||«i)M|W) 

»-£.  + («-£.))" 

' 4 (£t4£|-!1-£j  + (Ej-£|)f'(£i+£|-t. 
1 (il\\cd){kb\\al)(cJ\\kj) 

:-££  + (£.-£.))" 

2 (£.  + £,  - £r  - ££  + (£>  - £»))"’(£>  +£*-£, 
1 (kl\\ad){id\\cj)(cb\\kl) 

- £J  + (f,  - £.))* 

2 (£*  + £,  - £„  - ej  + (£j  — £6))A'(£i  + £(-£, 
{li\\da){kd\\d)(cb\\ki) 

1-£c+(£|-£«))A 

(e,  + £|  - f0  - £j  + (£,  - £*))"’(£>  + £*-£»- 

•£c + (£(-£.))" 

2 (£1  + £|  - £e  - ££  + (£j  - £i))A  (£j  + £J  - £c 

1 <il|M(M||«j)(crf||*l) 

- U + (f.  - c.)f 

2 (=,  + £,-£,.-££  + (£j  - £»))"'(£»  + £|  - ft 
1 

- £j  + (£,  - e.)f 

2 (£t  + £,-£„-££  + (£)  - ££»"'(££  +£(-£„ 
1 (iI|M)(d||j<>)(M|l*f} 

,-££+(£■-  £«))A 

(ft  + £(-£<-£,/  + (e>  - £i>)W'(f*  + £|  - £4  - £j  + (£i  - £o))f 


(kl\\ca)[ic\m(db\\lj) 


(£t4£|-£, 

1-£,  + (tr£/l£,  + t(-tt 
(i/||arf)(dH|je)(e6|HI) 

- ££  + (£■  “ £«))A 

(£!  + £.-£« 

, - ££  4 (£,  - £i))A"(£»  4 £t  - £i 

-£C  4 (£,-£.))" 

md«)m\d){cd\\kj) 

, - £j  + (e,  - «*))"'(«/  + =‘*  “ £' 

r — Id  4 (£i  - £.))* 

<fci||crf>(c/|K-a>(M||j/> 

{£,  4 £t  - £ 

, - ££  4 (£,-  - £*))A'(c*  4 £|  - £| 

|-££  4 (£.-£.))" 

(*/||a<0(ii||c/){a»||*j) 

o - ££  4 (£,  - Ci)f  (c,  4 £*  — ! 

:t-£J  4 (£,-£.))* 

(H||cd){cl|be)(W|U-0 

(£,4f*-i 

!,  - £j  4 (sj  - *»))"'(**  + £|  - ‘ 

:4-££4(£,-£.))A 

{/>||<fa)(»||Q-)W||M) 

U -££  4 (£j  -£»))*'(££  4 £|-< 

:t-££4(£,-£.))A 

, 

{H||cd)(ci||tn)(dHIO') 

v 

Dipole  Moment  Matrix  Elements 


Wl  = (X'klV')eV  + <x|ylv")e*  + (xl-'IWe. 
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<0|x|0)  = (*|x|t)  = xu 

(O|x|?)  = (i|x|0)  = x„ 

<®l*l«)  = 0 <289> 

<0|x|A'>  = 0 
(0|x|°i  A')  = 0 

(0|x|f)  = (i|x|o)  = x„  (290) 

(0|x|;f)  = 0 (291) 

(0|x|A')  = 0 (292) 

(0|x|“;  A')  = 0 (293) 

(‘|x|«)  = «i>4.i(0|x|0)  + 4i>(6|x|o)  - 6.t{i\x\j)  (294) 

($kl*i)  = Mj*(' |x|<0  + 4w4,i(t|x|c) 

(295) 

- - 4m4,j,(/|x|c) 
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(296) 


(297) 


(ifn  MS?)  = - <ml<n*)(0k|0) 

+ {tmk&nl  - iml^nlrt^/ctskW  + V(/Mc>+ 

-«tfM*W-<W/l*W]  (298) 

- (6/c< ,i  - + ft.lW*lm)+ 

-«„1,(i|*|n)-«.i(I|x|m)) 


_iuk — (,i,w  _ — ('-'»||mt>  ,, 

+ €,)'  ' Um  +£»-£/-  H Y 
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(/all"") 


' "(«..  + *-«/-« 

'g  4 (£i  — £«))”  " ' 

' S, . (fg\\mk) 

. Sim  (/Jilin) 

L (ffm  + £*-£/-  £,  + (£,  - ■ 

£«))*  (£4  4 £„-£/-£,  4 (£,  - £„))" 

(/c||mn) 

Sat  (cjllntn) 

l(£m4£„-£/-£c4(£,-. 

ta))K  (£„  4 £„-£<-£,  4 <£,  -£„))* 

[ (MM-) 

SimSat  (cjllin) 

!.■))*  («  + !.-£■-!»  + («-£.))'' 

f fcn*.,  (Mini) 

SmSat  (ej||im) 

l(£J  + £»-£/-£c  + (£i-£ 

.))"  (£r.  4 £4  -£.-£,  4 (£.  -£.))•'’ 

}«*M 

}(*l*|c) 

Wl«) 

(300) 


(A”|x|A’>  = (0|x|0> 


(301) 


rvr^T, 


(J;  AT'  |i|?;  A')  = (J;A"|f;  AT)(0|x|0)+ 


(303) 


(£m  + f,  - £C  - + (c>  - £*))'V'  ' Il|m)  (e*  + £1  - Ec  - U + <£.  - 

i <«||/d)  ,nik)  Ml**) ,, 

+ 2{ct  + c,-el-ei  + (c]-el))!<'  wl  ' (ft  + £/  - e, - ej  + (e;  - e.))A 


1^  ' 7 (ft  + £|  - £e  - f J + (£i  - £«))' 

_ u[xli)  MW 

l))-v  ’ 17  ' (ft  + £|  - £c  - £J  + (£.  - £«))■ 

,, W) 


"(££  + £. -£«-£/  + (£;  -ft))''' 


' (El  + e,  - £c  - ££  + (£|  - 
' (£t+£,-£c-££  + (£.-£- 


WM 


<‘w*> 

Ml«> 


4 (ft  + £|  - £c  - ej  + (Ej  - £i)); 

i mod)  ...  . . Wl«> 

2 (£|.  + £i-£.-£j + (£,-£»))''■■  ' ' (£»  + £!-£« - *4 + (£.■-£.))* 

' mi r m«i esa 

r(/W)7i 


(H||/«) 


+ 2 (£t  + £,-£/-£.  + (£j  - El))*'  |I|C'  (ft  + £|  - Et  - ft  + <£.  “ «.))V 

<*« F(%|,)S «- 


(ft  + £,  - £c  - £.  + (£,  - £&))'V" 


' (££ + £!-£<-££ + (£.  -£.»■' 
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i m\<*) 

5 (ft  4 e(  - £e  - £„  4 (e,-  - eh))K 

_1 (*'IM 

” 2 (£t  + £,  - £c  - ££  + (£,  - £»))* 

(fill*-) 

(£,  + £,-£t -£.  + (£, -£.))"' 
(fa'IW 

(£*  + £,  - £.  - £J  4 (£j  - £l));V" 


Ml*') 

£«-££+  (£.  ~ £« 
Ml*j) 

£e  - ££  + (£■  - £« 

<efc|l*J> 

£e  - ££  + (£i  - £.) 

Ml*j> 

:<-££  + (£. -£«): 


W 


l)-v 


(£*  + £i  — £< 
_(£t+£|-£, 

+ (£*  + £i-e« 


:-£.  + (£,-£.)r  V['UI 
fif! . 

:-£«  + (£>-  £*))V  ' ' ' ' I 


(#0 

(£4  4 £|  - £t  - £l  4 (£,  - £o)VV 

<cMI*i> _ 

(£l  4 £,  -£,-4  + (£.  - £.)),v 

(<*11  *;) 

(£t4£, -£c-££4(£,-£„))‘V 

Ml*;) , 

(e»  4£,-4-£d4  (£,  - £.))A' 


A great  simplification  can  be  achieved  if  spin  adaptation  is  included,  then  one 
would  obtain  a multi-electron  basis  set  for  singlets 

2.  If)  = S“'|0>;  s;'  = 4 «i.» 

3.  |JJ)  = SfSj"|0) 

4.  lA>  = (g,_lfl.)«Qt'l°)l  N=l,M 

5'  °-  = ir.-«.7Q^l‘?v|0):  A’  = I'M 
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r<?V|0);  N = l,M 


6.  Q = i - (|0>(0|  - ifXfl  - Iffxsfl) 


All  calculations  were  performed  using  1NDO/S  Model  Hamilionian  181,  82, 
101]  with  minimum  basis  sei. 

In  the  tables  17  lo  22  we  present  the  transition  energies  and  oscillator  strengths 
for  water,  formaldehyde,  and  acetylene.  We  compare  our  variational  results  (see 
table  16)  with  Cl-S.  Cl-SD.  CI-SDT  and  CI-F. 


NEXP1-ST_1  + doubles  (ij>>(a.b)  (|0),RV|0),  )|“)),{f?etiC?l'|0)],  {!#„))) 

aSfSgSb,  ifimSj,  i£n£j 

NEXP1-ST_0  + doubles  (ij)=>(a.b)  {|0), yjv'|0), {|f>}, {a*iflV|0)},{|4' )}} 

a’<f<g<b,  iSmSj.  i < n < j 

Clearly,  in  all  cases  Cl-SDT  gives  results  very  close  to  the  Cl-F  results. 
Cl-S  also  performs  quite  well  giving  the  right  order  of  the  states  except  for 
acetylene.  CI-SD  greatly  improves  the  transition  moments  and  yields  better 
oscillator  strengths. 


Table  16:  Description  of  the  NEXP  basis  sets  used. 


Method 

NEXPI-S 

NEXPl-ST.l 

NEXP1-ST_0 


Mullielcctron  Basis 


<|0),ftV'|0>,{|?>}} 

]|0),flvq0Mi?)},{«o'.v|0)}) 

)|0),flV'|0),{|f)},(a'iflV'|0)}) 


Ill 

NEXP1-S  presents  the  same  trends  as  CI-S,  except  with  larger  transition 
energies,  a consequence  of  the  fact  that  fiV|0)  correlates  the  ground  state.  This 
is  partly  corrected  by  the  inclusion  of  average  triple  excitations.  The  explicit 
inclusion  of  some  selected  double  excitations  yields,  in  general,  to  the  worst 
qualitative  results;  e.g.  we  note  an  intruder  slate  problem  for  NEXP1-ST_0  + 
doubles  (4,6)=>(7.8)  in  formaldehyde  and  an  inversion  of  stales  for  NEXP-ST_1 
+ doubles  (4,6)=>f7.8),  again  in  the  case  of  formaldehyde. 

For  acetylene,  NEXP1-ST_0  gives  the  same  ordering  of  the  slates  as  the 
Cl-SDT  for  the  first  three  transitions.  This  is  quite  remarkable  because  there  is 
a change  in  symmetry  from  II  to  E in  transitions  two.  three  and  four.  This  is 
not  destroyed  by  the  inclusion  of  double  excitations,  although  it  is  not  sufficient 
to  correct  the  order  of  the  fourth  transition  (— ).  The  NEXP1-ST_1  scheme  is 
related  closely  to  NEXP1-S  and  consequently  to  Cl-S.  NEXP1-S  presents  the 
same  ordering  as  Cl-S.  and  the  inclusion  of  triples  changes  the  nature  of  the 
founh  transition  from  Ilg  to  E„;  this  result  is  not  altered  by  the  inclusion  of  the 
selected  double  excited  states. 

For  acetylene,  the  inclusion  of  triples  results  in  belter  agreement  with  CI- 
SDT.  The  NEXP1-ST_0  scheme  scheme  also  performs  better  for  the  first  three 
transitions  of  formaldehyde.  This  suggests  that  the  single  excited  states  are  not 
very  sensitive  to  the  relative  weight  of  the  triple  excitations. 
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In  summary,  ihe  NEXP  spectroscopic  methods  presented  here  really  do  not 
perform  as  well  as  we  had  hoped,  although  for  most  systems  they  are  quite  good 
for  the  first  few  transitions.  Cl-S  and  CI-SDT  are  good  theories  for  the  lowest 
lying  states,  although  CI-SDT  becomes  prohibitively  expensive. 
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Table  20:  Excitation  Energies  and  Oscillator  Strengths  for  formaldehyde, 
dnc=  1,1020  A.  dco=l*2100  A.  «Hat=121.1°,  Eitf=-24.446548  Hanrees. 


Ground 

State 

(Hartrees) 

(Oscilator 

Strength) 

Excited 

State 

(Error) 

Transition 

Energy 

(cm1) 

(Oscilator 

Strength) 

Excited 

(Error) 

Energy 

(Oscilator 

Strength) 

State 

(Error) 

Transition 

Energy 

(Oscilator 

Strength) 

(Error) 

NEXP1-S 

-.096747 

41439.5 

(f.b.s) 

*2 

(127%) 

76248.9 

(.0241) 

b2 

(31%) 

96584.9 

(.2467) 

Is) 

(41%) 

118373.4 

(.0000) 

B| 

(50%) 

NEXP1- 

ST_0 

-.096983 

29812.1 

(f.b.s) 

Ay 

(64%) 

66523.1 

(.0216) 

By 

(14%) 

86962.8 

(.2710) 

(27%) 

102814.6 

(.0001) 

Bi 

(31%) 

NEXP1- 

ST_1 

-.097924 

30127.4 

(f.bj) 

A2 

(65%) 

66799.7 

(.0209) 

By 

(15%) 

86161.8 

(.2128) 

(26%) 

103533.4 

(.0002) 

Bt 

(32%) 

NEXP1- 

ST_0 

♦doubles 

(4,6)=>(7,8 

27556.0 

(f.b.s> 

(51%) 

64504.9 

(.0223) 

B2 

(11%) 

85309.5 

(.1595) 

(25%) 

86286.1 

(f.b.s) 

Ay 

NEXP1- 

ST.I+ 

doubles 

(4,6)=>(7,8 

12460.9 

(f.b.s) 

A2 

(32%) 

55384.8 

(.0201) 

By 

(5%) 

68186.4 

(.0041) 

(13%) 

80088.0 

(.0376) 

A, 

(17%) 

Cl-F 

-.107694 

18213.8 

A2 

58254.5 

By 

68240.8 

78592.5 

Bi 
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CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


In  chapter  3,  we  analyze  a Krylov-like  expansion  which,  although  it  can  be 
generalized,  we  use  to  expand  the  double-excited  states.  The  variational  energy 
obtained  by  using  the  basis  set  {|0),{flJsV|0)))  is  stable  at  all  intemuclcar 
distances  and  captures  over  95%  of  the  correlation  energy  obtained  by  CI-D.  This 
is  particularly  interesting  since  our  method  is  an  N5  method,  while  CI-D  is  N6, 
where  N is  the  number  of  basis  functions.  We  also  showed  the  connection  of  these 
functions  to  the  first  order  Brillouin-Wigner  Perturbation  Theory  wavefunclion. 
This  analysis  suggested  an  effective  non-iterative  method  to  obtain  the  Brillouin- 
Wigner  energies,  a method  which  converges  quickly  when  the  HOMO-LUMO 
gap  is  even  modest.  Our  non-iterative  method  of  solving  Brillouin-Wigner 
Perturbation  Theory  differs  from  that  of  Bartlett  and  Brandas  [19-23)  in  the  choice 
of  the  functions  used.  The  method  of  Bartlett  and  Brandas  functions  depends  on 
element  of  matrix  of  positive  powers  of  the  Hamiltonian,  while  our  functions 
depends  on  positive  powers  of  the  Rayleigh-Schroedinger  resolvent. 

In  chapter  4,  we  obtain  an  NJ.  method  for  polarizabilities  that  yields  good 
results  when  compared  with  CI-F.  This  method,  called  Modified  CI-S  Without 
Triples  (MCI-S_WT).  seems  particularly  suitable  to  applications  for  large  systems. 
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In  chapter  5,  we  present  two  ways  to  account  for  triple  excitations.  Both 
methods  seem  to  balance  the  average  inclusion  of  double  excitations,  resulting  in 
rather  good  transition  moments  (±10%  error  with  respect  to  CIF),  and  in  a good 
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